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Preface 


The purpose of this publication is to describe the 
FORTRAN library subprograms and their use in either 
a FORTRAN or an assembler language program. The 
body of the publication describes the mathematical 
subprograms (which perform computations) and the 
service subprograms (which perform testing and 
utility functions). This information is intended pri- 
marily for the FORTRAN programmer; Appendix E is 
intended for the assembler language programmer. 
Additional appendixes contain algorithms (the method 
by which a mathematical function is computed ), per- 
formance statistics, descriptions of interruption and 
error procedures, storage estimates, and sample stor- 
age printouts. 

The reader should be familiar with one of the follow- 
ing publications: 

IBM System/360 FORTRAN IV Language, Form 
C28-6515 

IBM System/360 Basic FORTRAN IV Language, 
Form C28-6629 

IBM System/360 Operating System: Assembler 
Language, Form C28-6514 

In addition, references are made within this publica- 
tion to information contained in the following publi- 
cations: 

IBM System/360 Principles of Operation, Form 
A22-6521 

IBM System/360 Operating System: Control Pro- 
gram Services, Form C28-6541 

Standard mathematical notation is used in this publi- 
cation. The reader is expected to be familiar with this 
notation and with common mathematical terminology. 


SECOND EpiT10n (November 1966 ) 


This publication is a major revision of Form C28-6596-0, and 
incorporates information released in Technical Newsletters 
N28-2151 and N28-2173. Significant changes have been made 
to support the Operating System FORTRAN IV Library in 
addition to the FORTRAN IV (E) Library. 

Significant changes or additions to the information contained 
in this publication will be reported in subsequent revisions or 
technical newsletters. 
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The Operating System/360 FrorTRAN tv library is com- 
prised of two types of relocatable subprograms: mathe- 
matical subprograms and service subprograms. The 
mathematical subprograms correspond to a subpro- 
gram defined by a FUNCTION statement in a FORTRAN 
source module. These subprograms always return one 
answer (function value) to the calling module. The 
service subprograms correspond to a subprogram de- 
fined by a SUBROUTINE statement in a FORTRAN source 
module. These subprograms may or may not return a 
value to the calling module. 

Calls to the library subprograms are either at the 
programmer's request or in response to program re- 


Introduction 


quirements. All calls are processed by the linkage 
editor, which takes the subprograms from the library. 
The library subprograms are then combined by the 
linkage editor with the calling module (either an 
object or a load module) into another load module 
which is ready for execution. 


The library subprograms may be called in either a 
FORTRAN or an assembler language program. The next 
two sections of this publication contain calling infor- 
mation for the FORTRAN programmer; Appendix E con- 
tains calling information for the assembler language 
programmer. 
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Mathematical Subprograms 


The mathematical subprograms supplied in the 
FORTRAN library perform computations frequently 
needed by the applications programmer. The mathe- 
matical subprograms are called in two ways: explicitly, 
when the programmer includes the appropriate entry 
name in a source language statement (see Table 1); 
and implicitly, when certain notation (e.g., raising a 
number to a power ) appears within a source language 
statement (see Table 6). 

The following text describes the individual mathe- 
matical subprograms and explains their use in a 
FORTRAN program. Detailed information about the 
actual method of computation used in each subpro- 
gram, the performance of the subprogram, interrup- 
tion and error procedures, and storage estimates can 
be found in the appendixes of this publication. 


Explicitly Called Subprograms 


Each explicitly called subprogram performs one or 
more mathematical functions. Each mathematical 
function is identified by a unique entry name that 
_ differs from the name of the subprogram. 

A subprogram is called whenever the appropriate 
entry name is included in a FoRTRAN arithmetic expres- 
sion. The programmer must also supply one or more 
arguments. These arguments follow the entry name 
and are separated by commas; the list of arguments is 
enclosed in parentheses. 

For example, the source statement 


RESULT = SIN (RADIAN) 


causes the IHCsscN subprogram to be called. The sine 
of the value in RADIAN is computed and the function 
value is stored in RESULT. 

In the following example, the mHcsserT subprogram 
is called to compute the square root of the value in 
AMNT. The function value is then added to the value 
in STOCK and the result is stored in ANs. 


. ANS = STOCK + SQRT (AMNT) 


The explicitly called subprograms are described in 
the tables which make up the rest of this section. 
These tables show the general function, subprogram 
name, the FORTRAN library that contains the subpro- 
gram, definition, entry name(s), argument information, 
type of function value returned, and assembler re- 
quirements. The following column headings are used 
in the tables: 


6 


General Function: This column states the nature of 
the computation performed by the subprogram. 

Subprogram Name: This column gives the module 
name of the subprogram. (The ForTRAN library is 
assigned an internal code of mc which precedes each 
subprogram name. ) 

Subset: This column indicates those subprograms 
that belong to the ForTRAN Iv (£) library. Unless 
otherwise indicated, all subprograms that belong to 
the £ library also belong to the rorTRAN Iv library. 

Definition: This column gives a mathematical equa- 
tion which represents the computation. An alternate 
equation is given in those cases where there is another 
way of representing the computation in mathematical 
notation. (For example, the square root can be repre- 
sented either as y= x or y = x!/2,) The definition 
for those subprograms that accept complex arguments 
contains the notation z = x,+ Xgi. 

Entry Name: This column gives the entry name that 
the programmer must use to call the subprogram. A 
subprogram may have more than one entry name; the 
particular entry name used depends upon the compu- 
tation to be performed. For example, the micsscn sub- 
program has two entry names: cos and sin. If the 
cosine is to be computed, entry name cos is used; if 
the sine is to be computed, entry name sin is used. 

Argument Number: This column gives the number 
of arguments that the programmer must supply. 

Argument Type: This column describes the mode 
and length of the argument. INTEGER, REAL, and 
COMPLEX represent the type of number; the notation 
*4, *8, and *16 represent the size of the argument in 
storage locations. 


Note: In ForTRAN Iv (E), a real argument corre- 
sponds to the REAL *4 argument, and a double- 
precision argument corresponds to the REAL *8 argu- 
ment. Complex arguments cannot be used with a 
FORTRAN Iv (E) compiler. 


Argument Range: This column gives the valid range 
for an argument. If the argument is not within this 
rane, an error message is issued and execution of this 
load module is terminated. Appendix C contains a 
description of the error messages. 

Function Value Type: This column describes the 
type of function value returned by the subprogram. 
The notation used is the same as that used for the 
argument type. 


Table 1. Explicitly Called Mathematical Subprograms 


















Logarithmic and exponential 
subprograms (described in Table 2) 


Exponential 


Square root 





Trigonometric subprograms Arcsine and arccosine 


(described in Table 3) 


Arctangent 


Sine and cosine 








Hyperbolic function subprograms 
(described in Table 4) 








Miscellaneous subprograms 
(described in Table 5) 


Error function 





Modular arithmetic 


Truncation 








Assembler Requirements: This column gives the reg- 
isters used by the subprogram and the minimum save 
area that the assembler language programmer must 
supply. For example, the assembler requirements for 
the rHccssQeT subprogram are: 

registers 0, 2(4) 
save area 9F 
This information specifies that: 


1. The function value is found in floating-point reg- 


isters 0 and 2. 


Tangent and cotangent 










Gamma and log-gamma 


Maximum and minimum value 





General Function Specific Function Subprogram Name Entry Name(s). 


Common and natural logarithm 


IHCCLLOG* CDLOG 





IHCCSLOG* CLOG 

IHCLLOG DLOG, DLOG10 
IHCSLOG ALOG, ALOG10 
IHCCLEXP* 

IHCCSEXP* 

IHCLEXP 

IHCSEXP 

IHCCLSQT* CDSQRT 
IHCCSSQT* CSQRT 
IHCLSQRT DSQRT 
IHCSSQRT SQRT 
IHCLASCN* DARSIN, DARCOS 
IHCSASCN* ARSIN, ARCOS : 
IHCLATAN DATAN 
IHCLATN2* DATAN, DATAN2 
IHCSATAN ATAN 
IHCSATN2* ATAN, ATAN2 
IHCCLSCN* CDSIN, CDCOS 
IHCCSSCN* CSIN, CCOS 
IHCLSCN DSIN, DCOS 
IHCSSCN SIN, COS 
IHCLTNCT* DTAN, DCOTAN 
IHCSTNCT* TAN, COTAN 


Hyperbolic sine and cosine IHCLSCNH* DSINH, DCOSH 
IHCSSCNH* SINH, COSH 

Hyperbolic tangent IHCLTANH | DTANH 
IHCSTANH TANH 

Absolute value IHCCLABS* 
THCCSABS* 





CDABS 

CABS 
IHCLERF* DERF, DERFC 
IHCSERF* ERF, ERFC _ 
IHCLGAMA* DGAMMA, DLGAMA 
IHCSGAMA* GAMMA, ALGAMA 









IHCFMAXD DMAX]1, DMIN1 
IHCFMAXI AMAXO, AMINO, MAXO, MINO 
IHCFMAXR AMAX]1, AMINI, MAX], MIN1 


IHCFMODI MOD 
IHCFMODR AMOD, DMOD 
IHCFAINT AINT 
IHCFIFIX IDINT, INT 


2. Floating-point register 4 is used for intermediate 
computation. 
3. The save area must be at least nine full-words in 
length. 
Detailed information for the assembler language pro- 
grammer is given in Appendix E. 


Note: In the following tables, the approximate value 
of 218 « x is .82354966406249996D + 06; the approximate 
value of 25° « = is .35371188737802239D + 16. 
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Table 2. Logarithmic and Exponential Subprograms 


























































A : a ae ’ eo eee. Assembler 
General Subprogram | Sub- Entry ESUIMEDS Value Require- 
Function Name set Definition Name Type’ | —sRange Type’ ments 







registers 
0, 2 
save area 


8F 


registers 
0, 2 

save area 
8F 


z 40 + Oi 
See Note 3 


CDLOG 






y=Ppv log. (z) 
See Note 2 


Common IHCCLLOG 
and natural 


logarithm 
















complex 
*16 


complex 

*16 
z~0+ 0i complex 
See Note 3 *8 


*8 









CLOG 


° 


y=pPpv log. (z) 
See Note 2 


IHCCSLOG 


| ae 
° 


IHCLLOG Yes y=logex or DLOG real *8 x >0 real *8 registers 
=Inx 0 (2) 
save area 
OF 


registers 
0 (2) 
save area 
OF 


DLOGI10 


IHCSLOG Yes y=logex or ALOG registers 
y=lInx 0 (2) 
save area 


ia - 
~ —_ iin 


5F. 


registers 
0 (2) 
save area 


5F 


registers 
0, 2 

save area 
8F 
registers 
0, 2 

save area 


8F 


registers 


0 (2) 


Exponential | IHCCLEXP | No 


CDEXP complex | x: S 174.673 complex 
*16 Ixe] < (2%) | *16 
complex | x SS 174.673 complex 

*8 \xo| << (2° ar) | *8 

real *8 x S 174.673 real *8 
save area 
oF 

real *4 x S 174.673 real *4 register 0 
save area 
12F 

1 complex | any complex complex | registers 

*16 argument *16 0, 2 (4) 
See Note 3 save area 
OF 
complex | any complex complex 
*8 argument *8 
See Note 3 
Sea 
Viel ll 
NOTEs: 


registers 

0, 2 (4) 

save area 
OF 
registers 

1. In FORTRAN IV (E), a real argument corresponds to the REAL *4 argument, and a double-precision argument corresponds to 

| the REAL *8 argument. 

2. PV = principal value. The answer given is from that point where the imaginary part (y2) lies between —7 and +7. More specifi- 

cally: —w < ya S 1, unless x1 < 0 and x2 = —O, in which case, ye = —7. 

3. Floating-point overflow can occur. 





THCCSEXP CEXP 





IHCLEXP Yes DEXP 






THCSEXP Yes 


Square root | IHCCLSQT | N CDSQRT 


No CSQRT 


IHCCSSQT 








IHCLSQRT | Yes DSQRT 
0 (2, 4) 
save area 


5F 


registers 
0 (4) 

save area 
5F 









IHCSSQRT Yes SQRT 





| ae 
° 









Table 3. Trigonometric Subprograms 


Subprogram | Sub- Entry 
Name set Definition Name 


| THCLASCN | No y = arcsin (x) DARSIN 
y = arccos (x) DARCOS 


y = arcsin (x) ARSIN 


IHCSASCN | No 
_ 













General 
Function 


rm Function | Assembler 
rgument(s) Value Require- 
No. Type’ Range Type’ ments 

















Arcsine and 1 | real *8 real registers 
arccosine *8 (in 0 (2, 4) 
radians) | save area 
13F 













real 
*8 (in 
radians ) 


real *8 






registers 
0 (2, 4) 
save area 


13F 


















registers 
0 (2, 4) 
save area 
10F 


real 
*4 (in 
radians ) 


real *4 


















registers 
0 (2, 4) 
save area 
10F 


real 
*4 (in 
radians ) 


real *4 

























Arctangent | IHCLATAN real *8 any real real registers 
argument *8 (in 0 (2, 4, 6) 
radians) | save area 


5F 


_ 



















IHCLATN2 y = arctan (x DATAN real *8 any real real registers 
argument *8 (in 0 (2, 4, 6) 
radians ) | .save area 
5F 

















" Xi DATAN2 real *8 any real real registers 
y = arctan | 7° arguments *§8 (in 0 (2, 4, 6) 
(except 0, 0) radians) | save area 
5F 


























IHCSATAN real *4 any real real registers 
argument *4 (in 0 (2, 4, 6) 
radians) | save area 


5F 


- 


IHCSATN2. y = arctan (x) ATAN 


















registers 
0 (2, 4, 6) 
save area 


oF 


real 
*4 (in 
radians ) 


real *4 any real 


argument 


























ATAN2 real *4 any real real registers 
arguments *4 (in 0 (2, 4, 6) 
(except 0, 0) radians) | save area 





5F 


IHCCLSCN 






CDSIN 


— 


complex 
*16 (in 
radians ) 


registers 
0, 2 (4) 
save area 


OF 


[x1]< (2° * 3) 
Ixs| S 174.673 





Sine and 
cosine 























complex 
*16 









CDCOS complex | |xi|<(2 + 7) complex | registers 
“16 (in | |xo| 174.673 | *16 0,2 (4) 
radians ) save area 

OF 


IHCCSSCN y = sin (z) CSIN 





|x1|<( 2° * ar) 
xe| = 174.673 






complex 
*8 (in 
radians ) 





complex | registers 
“8 0, 2 (4) 
save area 
OF 
complex | |x:|<(2% +7) complex | registers 
*8 (in lxe] 174.673 | *8 0, 2 (4) 
radians ) save area 
OF 
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y = cos (z) CCOS 





Table 3. Trigonometric Subprograms (Continued ) 


Subprogram 
Name 


General 


Function Definition 
Sine and IHCLSCN 
cosine 


(continued ) 


y = sin (x) 


IHCSSCN 
Tangent and| IHCLTNCT | No 
cotangent 

IHCSTNCT 


y = tan (x) 


Notes: 


y = cotan (x) 


y = cotan (x) 


DCOTAN 


es 


No. 


Argument(s) 


real 
*8 (in 
radians ) 


real 
*8 (in 
radians ) 


real 
*4 (in 
radians ) 


real 
*4 (in 
radians ) 


real 
*8 (in 
radians ) 


real 
*8 (in 
radians ) 


real 
*4 (in 
radians ) 


real 
*4 (in 
radians ) 


|x]<(2% © ar) 


Ix} <(2 +) 


|x]<(2 © 3) 
See Note 3 


|x|<( 2° « 3) 
See Note 3 


[x|<<( 2" * 3) 
See Note 3 


xi 27 ar) 
See Note 3 


Function 
Value 
Type’ 


Assembler 
Require- 
ments 


registers 
0 (2, 4) 
Save area 
5F 


registers 
0 (2, 4) 
save area 
5F 


registers 
0 (2, 4) 
save area 
5F 


registers 
0 (2, 4) 
save area 


5F 


registers 
0 (2, 4, 6) 
save area 
5F 


registers 
0 (2, 4, 6) 
save area 
5F 


registers 
0 (2, 4) 
save area 


5F 


registers 
0 (2, 4) 
Save area 


5F 


. In FORTRAN IV (E), a real argument corresponds to the REAL *4 argument, and a double-precision argument corresponds to 


the REAL *8 argument. 


2. Instead of the IHCLATAN and IHCSATAN subprograms contained in the FORTRAN IV (E) library, the FORTRAN IV library 


contains the IHCLATN2 and IHCSATN2 subprograms. 


. The argument for the cotangent functions may not be near a multiple of +; the argument for the tangent functions may not be near 


an odd multiple of 1/2. 
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Table 4. Hyperbolic Function Subprograms 






















A t(s) Function | Assembler 
General Subprogram | Sub- Entry TEUMENECS Value Require- 
Function Name set Definition Name No. Type* Range Type’ ments 





DSINH 1 real *8 registers 
0 (2, 4) 
save area 


OF 


THCLSCNH 


IHCSSCNH 





Hyperbolic 
sine and 
cosine 












_ 
Ix|<174.673 
real *4 | |x| <174.673 real *4 
save area 
OF 
real *4. | |x|<174.673 real*4 ‘| registers 
0 (2, 4) 
save area 
OF 
any real real *8 | registers 
argument 0 (2) 
save area 
5F 
any real real *4 registers 
argument 0 (2) 
save area 
5F 
NoTE: 


1, In FORTRAN IV (BE), a real argument corresponds to the REAL *4 argument, and a double-precision argument corresponds 
to the REAL *8 argument. 


DCOSH real *8 





registers 
0 (2, 4) 
save area 
OF 









registers 
0 (2, 4) 


SINH 







COSH 








IHCLTANH DTANH real *8 







Hyperbolic 
tangent 





IHCSTANH TANH real *4 














Table 5. Miscellaneous Mathematical Subprograms 


Sub- Entry Argument(s ) 
set Definition Name Type* 


y=(|z| = (x2+x.")? 1 | complex 
*16 

















Assembler 
Require- 
ments 


Function 
Value 
Type’ 















Subprogram 
Name 


General 
Function 



























Absolute IHCCLABS any complex registers 
value argument 0, 2 (4) 
See Note 5 save area 






8F 






























‘IHCCSABS any complex real *4 registers 
argument 0, 2 (4) 
See Note 5 save area 
8F 


THCLERF 






real *8 
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registers 
0 (2, 4, 6) 
Save area 
11F 





any real 
argument 


Error 
function 



































registers 
0 (2, 4, 6) 
save area 

11F 


any real 
argument 






ae 
even es 
y=1l-—erf (x) 






aa 
*8 


Table 5. Miscellaneous Mathematical Subprograms (Continued ) 







































General Subprogram | Sub- Entry TeHINERtS Value Require- 
Function Name set Definition Name No. Type’ Type’ ments 



































Error IHCSERF real *4 any real real *4 registers 
function argument 0 (2, 4, 6) 
(continued ) save area 
l1F 
ERFC real *4 any real real *4 registers 
argument 0 (2, 4, 6) 
save area 
11F 


Gamma and DGAMMA real *8 x > 2° and real *8 registers 

log-gamma x < 57.5744 0 (2, 4, 6) 
save area 
11F 


IHCLGAMA 


THCSGAMA 
IHCFMAXD 


IHCFMAXI 


DLGAMA real *8 x > O and real *8 registers 
x << 4,2913 + 10" 0 (2, 4, 6) 
save area 
11F 





% GAMMA real *4 x > 2°* and real *4 registers 
y= fom e' du x <. 57.5744 0 (2, 4, 6) 
0 save area 
11F 


y=log. T (x) or ALGAMA real *4 x > Oand real *4 registers 
0 x <4,2913* 10" 0 (2, 4, 6) 

y=log. f u**e™" du save area 
0 11F 


IV i e ee i —_ 
bo 





Maximum real *8 
and 
minimum 
values 


DMAX1 register 0 
save area 


oF 


y=max (x1,... 





any real real *8 
arguments 
any real real *8 
arguments 





DMINI1 real *8 


IIV 
bo 


register 0 
save area 
OF 


y=min (x,..., Xn) 


y=log. T (x) or 
oe 

y=log. f u*te" du 
0 


No 
No 
No 
Yes 

es 





y=max (xi1,...,Xn) | AMAXO =2 | integer any integer real *4 register 0 
arguments Save area 
oF 
MAX0 =2 | integer any integer integer register 
arguments See Note 2 
save area 
OF 





AMINO 


IV 
bo 


register 0 
save area 


y=min (x1,..., Xn) integer 





Bs) 
3) 





any integer real *4 
arguments 










OF 
=2 | integer any integer integer register 
arguments See Note 2 
save area 
OF 
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Table 5. Miscellaneous Mathematical Subprograms (Continued ) 


Entry 
Name No. 









Argument(s) 
ype" 







Subprogram | Sub- 
Name set 


General 
Function 


we 


Definition Range 


































Maximurn IHCFMAXR y=max (x,...,Xn) | AMAXI =2] real *4 any real 

and arguments 

Minimum 

Values 

(continued ) MAX1 =>2] real *4 | any real 
AMINI 


arguments 


real *4 


IIV 
bo 


+5 Xu) 





y=min (x,.. any real 


arguments 


IV 
bo 


real *4 any real 


arguments 








Modular JHCFMODI MOD 


arithmetic 


y = x: (modulo x2) 
See Note 3 


X2 54 0 
See Note 6 


integer 


IHCFMODR y = x: (modulo xz) AMOD real *4 Xe ~ 0 
See Note 3 See Note 6 
y = x: (modulo xz) real *8 Xe 34 0 
See Note 3 See Note 6 


DMOD 
AINT 





IACFAINT y = (signx)¢n real *4 
where n is the largest 


integer < |x| 


Truncation any real 


argument 






IHCFIFIX See 
Note 


y = (signx) +n real *8 
where n is the largest 


integer < |x| 


any real 
argument 


real *4 


any real 
argument 


Notes: 


to the REAL *8 argument. 


2. The result is stored in general register 0. 


5. Floating-point overflow can occur. 


tion occurs. (A detailed description of the interruption procedure is given in Appendix C.) 














Function 


Value 
Type’ 


teger whose magnitude does not exceed the magnitude of *! is used. The sign of the integer is the same as the sign of i 
x 





real *4 
real *8 ; 


Assembler 
Require- 
ments 


register 0 
save area 


OF 


register 
See Note 2 
save area 


OF 





register 0 
save area 
OF 


register 
See Note 2 
save area 
OF 


register 
See Note 2 
save area 
OF 


register 0 
save area 
OF 


register 0 
save area 


OF 


register 0 
save area 


OF 


register 0 
save area 


OF 





register 
See Note 2 
save area 
OF 


1. In FORTRAN IV (E), a real argument corresponds to the REAL *4 argument, and a double-precision argument corresponds 


3. The expression x: (modulo xz) is defined as x: — [=] * x2, where the brackets indicate that an integer is used. The largest in- 
Xe 


Xa 


4, The coding that performs this function is out-of-line in FORTRAN IV (E) and in-line in FORTRAN IV. Out-of-line coding is 
taken from the FORTRAN library by the linkage editor and processed with the calling module. In-line coding is inserted by the 
FORTRAN compiler at the point in the source module where the function is referenced. This means that the in-line functions 
are available in FORTRAN IV by using the appropriate entry name but that they are not part of the library. 


6. If x2 = 0, then the modulus function is mathematically undefined. In addition, a divide exception is recognized and an interrup- 
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Implicitly Called Subprograms 


The implicitly called subprograms perform operations 
required by the appearance of certain notation in a 
FORTRAN source statement. When a number is to be 
raised to a power or when multiplication and division 
of complex numbers are to be performed, the FORTRAN 
compiler generates the instructions necessary to call 
the appropriate subprogram. For example, if the fol- 
lowing source statement appears in a source module, 


ANS = BASE**EXPON 
where BASE and EXPON are values of the form REAL *4, 


the IHCFRxPR subprogram is called by the FORTRAN 
compiler. 


The implicitly called subprograms in the FORTRAN 
library are described in Table 6. This table shows the 


Table 6. Implicitly Called Mathematical Subprograms 









Implicit 









Raise an integer 
to an integral 
power 


IHCFIXPI Yes 


Raise a real THCFRXPI Yes 
number to an 


integral power 


THCFDXPI 
Raise a real JIHCFPXPR 
number to a 
real power IHCFDXPD 


Raise a complex 
number to an 
integral power 


IHCFCDXI 


NOTEs: 


must be included. 


the REAL *8 argument. 
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General Subprogram | Sub- Function Entry’ Argument(s) Value Assembler 
Function Name set Reference* Name Type® Requirements 
Multiply and IHCCLAS No y = 21" Ze CDMPY# 2 | complex *16 | complex *16 | registers 0, 2 (4, 6) 
divide complex save area 5F 
numbers CDDVD# complex *16 | complex *16 | registers 0, 2 (4, 6) 
save area 5F 
IHCCSAS No CMPY# complex *8 complex *8 registers 0, 2 (4, 6) 
save area 5F 
CDVD# complex *8 complex *8 registers 0, 2 (4, 6) 
save area 5F 


FIXPI# 
FRXPI#. | 2 
FDXPI# 
FRXPR# 
FDXPD# 


FCDXI# 2 
IHCFCXPI No y=z**j FCXPI# 2 


1. This is only a representation of a FORTRAN statement; it is not the only way the subprogram may be called. 
2. This name must be used in an assembler language program to call the subprogram; the character # is a part of the name and 


3. In FORTRAN IV (E), a real argument corresponds to the REAL *4 argument and a double precision argument corresponds to 


general function, subprogram name, the FORTRAN 
library that contains the subprogram, implicit function 
reference, entry name, argument information, type of 
function value returned, and assembler requirements. 
The column headed “Implicit Function Reference” 
gives a sample source statement that might appear in 
a FORTRAN source module and cause the subprogram 
to be called. The rest of the column headings in Table 
6 have the same meaning as those used with the ex- 
plicitly called subprograms. ; 

The action taken within the subprogram depends 
upon the type of base and exponent used. Tables 7-10 
show the result of an exponentiation performed with 
the different combinations and values of base and ex- 
ponent. In these tables, I and J are integers; A and B 
are real numbers; C is a complex number. 











Function 

























i= integer 
j = integer 


register 0 
save area 18F 






real *4 register 0 
save area 18F 

real *8 register 0 
save area 18F 

real *4 register 0 
save area 18F 


register 0 
save area 18F 







a = real *4 
j = integer 


a = real *8 


bo 
op] 
i il 
et 
oo 
a 
* & 
-~ 


Ne] 








register 0 
save area 18F 










complex *8 register 0 


save area 18F 


Z = complex 
*8 
j = integer 






Table 7. Exponentiation With Integer Base and Exponent 


















Exponent (J) 







Base (1) Ts 6 


Compute the 
function value 

















Function value = 1 
ifI=1 

Otherwise, function 

value = 0 


Function 
value = 1 











Function 
value = 0 

















Function value= — 1 
ifI = —l and 
if J is an odd 
number 
Function value = 1 
if I = —1 and 
if J is an even 
number 
Otherwise, function 
value = 0 


Function 
value = 1 


Compute the 
function value 







Table 8. Exponentiation with Real Base and Integer Exponent 














Exponent (J) 


Function 


Compute the Compute the 
function value 


function value | value = 1 

























Function Error message Error message 
value = 0 THC242I or THC242I or 
THC2431 THC2431 
Compute the | Function Compute the 
function value value =] function value 


Table 9. Exponentiation with Real Base and Exponent 
Exponent (B) 
pase (A) Eom (BY) 
“oT p>o [| B-0 |  B<o | 
Compute the | Function Compute the 
function value| value = 1 function value 


Function Error message Error message 
value = 0 IHC244I or IHC244I or 
IHC2451 THC2451 










A> 








0 
0 



















A<0 | Error message | Function Error message 
THC253I or value = 1 THC253I1 or 
IHC2631 IHC2631 


Table 10. Exponentiation with Complex Base and Integer 
Exponent 


Base (6) 
C=R+Ri 












R > 0 and Compute the | Function Compute the 
Ri > 0 function value | value=1 + Oi| function value 
R > 0 and i 

Ri = 0 







Compute the | Function Compute the 
function value | value=1 + Oi] function value 
R > 0 and Compute the | Function Compute the 

Ri <0 function value | value=1 + Oi] function value 


R = O and Compute the | Function Compute the 
Ri > 0 function value | value=1 + Oi] function value 






R = O and 
Ri = 0 


Function value] Error message | Error message 
IHC246I or THC246I or 
IHC2471 IHC2471 


Compute the | Function Compute the 
function value | value=1 + Oi] function value 
Compute the | Function Compute the 
function value | value=1 + Oi} function value 
Compute the | Function Compute the 
function value | value=1 + Oi} function value 




















R < Oand 


Ri = 0 
R < Oand Compute the | Function Compute the 
Ri <0 function value | value=1 + Oi | function value 
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Service Subprograms 


The service subprograms supplied in the FORTRAN 
library are divided into two groups: one group tests 
machine indicators and the other group performs 
utility functions. Service subprograms are called by 
using the appropriate entry name in a FORTRAN source 
language CALL statement. 


Machine Indicator Test Subprograms 


The machine indicator subprograms (IHCFSLIT, 
IHCFOVER, and IHCFDVCH) test the status of pseudo in- 
dicators and may return a value to the calling program. 
When the indicator is zero, it is off; when the indica- 
tor is other than zero, it is on. In the following descrip- 
tions of the subprograms, i represents an integer 
expression and j represents an integer variable. 


IHCFSLIT Subprogram 


The 1HcFsLir subprogram is used to alter, test, and/or 
record the status of pseudo sense lights. Either of two 
entry names (SLITE or SLITET) is used to call the sub- 
program. The particular entry name used in the CALL 
statement depends upon the operation to be per- 
formed. 

If the four sense lights are to be turned oFF or one 
sense light is to be turned on, entry name SLITE is used. 
The source language statement is: 


CALL SLITE(i) 


where i has a value of 0, 1, 2, 3, or 4. 

If the value of i is 0, the four sense lights are turned 
off; if the value of i is 1, 2, 3, or 4, the corresponding 
sense light is turned on. If the value of i is not 0, 1, 2, 
3, or 4, an error message is issued and execution of 
this load module is terminated. (This error message 
is explained in Appendix C.) 

If a sense light is to be tested and its status recorded, 
entry name sLITET is used. The source language state- 
ment is: 


CALL SLITET (i,j) 
where: 


i has a value of 1, 2, 3, or 4, and indicates which 
sense light to test. 

j is set to 1 if the sense light was on; or to 2 if the 
sense light was off. 


If the value of i is not I, 2, 3, or 4, an error message 
is issued and execution of this load module is ter- 
minated. (This error message is explained in Appen- 
dix C.) 
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IHCFOVER Subprogram 
The IHCFOVER subprogram tests for an:exponent over- 
flow or underflow exception and returns a value that 
indicates the existing condition. After testing, the over- 
flow indicator is turned off. This subprogram is called 
by using the entry name OVERFL in a CALL statement. 
The source language statement is: 
CALL OVERFL (/) 
where: 
j is set to 1 if a floating-point overflow condition 
exists; to 2 if no overflow or underflow condition 
exists; or to 3 if a floating-point underflow condi- 
tion exists. A detailed description of each exception 
is given in Appendix C. 


IHCFDVCH Subprogram 


The rHcrpvcH subprogram tests for a divide-check 
exception and returns a value that indicates the exist- 


_ing condition. After testing, the divide-check indicator 
‘is turned off. This subprogram is called by using entry 


name DVCHK in a CALL statement. The source language 
statement is: 
CALL DVCHK (/) 
where: 
j is set to 1 if the divide-check indicator was on; or 
to 2 if the indicator was off. A detailed description of 
the divide-check exception is given in Appendix C. 


Utility Subprograms 

The utility subprograms perform two operations for 
the FORTRAN programmer: they either terminate execu- 
tion (IHCFExIT) or dump a specified area of storage 
(IHCFDUMP ). ' 


IHCFEXIT Subprogram 

The 1HCFExIT subprogram terminates execution of this 
load module and returns control to the operating sys- 
tem. (This subprogram performs a function similar to 
that performed by the stop statement.) The IHCFEXIT 
subprogram is called by using the entry name ExIT 
in a CALL statement. The source language statement is: 


CALL EXIT 


IHCFDUMP Subprogram 

The 1HcrpuMP subprogram dumps a specified area of 
storage. Either of two entry names (DUMP or PDUMP) 
can be used to call the subprogram. The entry name 


is followed by the limits of the area to be dumped and 
the format specification. The entry name used in the 
CALL statement depends upon the nature of the dump 
to be taken. 

If execution of this load module is to be terminated 
after the dump is taken, entry name pump is used. The 
source language statement is: 


CALL DUMP (a:, bi, fi, ..., dn, bn, fn) 

where: 
a and b are variables that indicate the limits of stor- 
age to be dumped (either a or b may represent the 
upper or lower limits of storage). 
f indicates the dump format and may be one of the 
integers given in Table 11. The formats available 
depend upon the compiler in use. A sample printout 
for each format is given in Appendix F. 


Table 11. The IHCFDUMP Subprogram Format Specifications 


FORTRAN IV (E) FORTRAN IV 


0 specifies hexadecimal 0 specifies hexademical 
4 specifies INTEGER 1 specifies LOGICAL *1 
5 specifies REAL 2 specifies LOGICAL *4 
6 specifies DOUBLE 3 specifies INTEGER *2 
PRECISION 4 specifies INTEGER *4 
5 specifies REAL *4 
6 

7 

8 

9 

























specifies REAL *8 
specifies COMPLEX *8 
specifies COMPLEX *16 
specifies literal 








If execution of this load module is to be resumed 
after the dump is taken, entry name ppump is used. 
The source language statement is: 


CALL PDUMP (a, bi, fh eee 5 Any ba, fn) 


where a, b, and f have the same meaning as explained 
previously. 


Programming Considerations 


A load module may occupy a different area of storage 
each time it is executed. To ensure that the appropri- 
ate areas of storage are dumped, the following conven- 
tions should be observed. 


Note: In the following examples, a is a variable in 
COMMON, B is a real number, and the array TABLE is 
dimensioned as: 


DIMENSION TABLE (20) 
If an array and a variable are to be dumped at the 


same time, a separate set of arguments should be 
used for the array and for the variable. The specifica- 


tion of limits for the array should be from the first 
element in the array to the last element. For example, 
the following call to the m1crpump subprogram could 
be used to dump TABLE and B in hexadecimal format 
and terminate execution after the dump is taken: 

CALL DUMP (TABLE (1), TABLE (20), 0, B, B, 0) 

If an area of storage in COMMON is to be dumped 
at the same time as an area of storage not in COMMON, 
the arguments for the area in COMMON should be given 
separately. For example, the following call to the 
IHCFDUMP subprogram could be used to dump the 
variables A and B in REAL *8 format without terminat- 
ing execution: 

CALL PDUMP (A,A,6,B,B,6 ) 

If variables not in comMMON are to be dumped, each 
variable must be listed separately in the argument 
list. For example, if rn, Pp, and Q are defined implicitly 
in the program, the statement 


CALL PDUMP (RB,R,5,P,P,5,Q,Q,5) 


should be used to dump the three variables. If the 
statement 

CALL PDUMP (R,Q,5) 
is used, all main storage between rR and Q is dumped, 
which may or may not include p, and may include 
other variables. 

If an array and a variable are passed as arguments 
to a subroutine, the arguments in the call to the 
IHCFDUMP subprogram in the subroutine should specify 
the parameters used in the definition of the subroutine. 
For example, if the subroutine sunt is defined as: 

SUBROUTINE SUBI (X, Y) 
DIMENSION X(10) 
and the call to sust within the source module is: 
DIMENSION A(10) 


CALL SUBI (A, B) 


then the following statement in the subroutine should 
be used to dump the variables in hexadecimal format 
without terminating execution: 


CALL PDUMP (X(1), X(10), 0, Y, Y, 0) 
If the statement 
CALL PDUMP (X(1), Y, 0) 


is used, all storage between A(1) and y is dumped, 
due to the method of transmitting arguments. 
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Appendix A. Algorithms 


Appendix A contains information about the computations used in the explicitly 
called mathematical subprograms. This information is arranged in alphabetical 
order, according to the module name of the subprogram. The entry names asso- 
ciated with each subprogram are given in parentheses after the module name. 


The information for each subprogram is divided into two parts. The first part 
describes the algorithm used; the second part describes the effect of an argument 
error upon the accuracy of the answer returned. 


The presentation of each algorithm is divided into its major computational steps; 
the formulas necessary for each step are supplied. Some of the formulas are 
widely known; those that are not so widely known are derived from more common 
formulas. The process leading from the common formula to the computational 
formula is sketched in enough detail so that the derivation can be reconstructed 
by any one who has an understanding of higher mathematics and access to the 
common texts on numerical analysis.! 


The accuracy of an answer produced by these algorithms is influenced by two 
factors: the performance of the subprogram (see Appendix B) and the accuracy 
of the argument. The effect of an argument error upon the accuracy of an answer 
depends solely upon the mathematical function involved and not upon the partic- 
ular coding used in the subprogram. 


A guide to the propagational effect of argument errors is provided because 
argument errors always influence the accuracy of answers whether the errors are 
accumulated prior to use of the subprogram or introduced by newly converted 
data. This guide (expressed as a simple formula where possible) is intended to 
assist users in assessing the effect of an argument error. 


The following symbols are used in this appendix to describe the effect of an 
argument error upon the accuracy of the answer: 





SYMBOL EXPLANATION 
q(x) The result given by the subprogram. 
f(x) The correct result. 
f(x) — g(x) The relative error of the result 
: pe given by the subprogram. 
8 The relative error of the argument. 
E | f(x) — g(x) | The absolute error of the result 
given by the subprogram. 
A The absolute error of the argument. 


The notation used for the continued fractions complies with the specifications 
set by the National Bureau of Standards. For more information, see Milton 
Abramowitz and Irene A. Stegun (editors), Handbook of Mathematical Functions, 
Applied Mathematics Series-55 (National Bureau of Standards, Washington, 
D.C.), 1965. 


1Any of the common numerical analysis texts may be used as a reference. One such text is F. B. Hildebrand’s 

Introduction to Numerical Analysis (McGraw-Hill Book Company, Inc., New York, N. Y., 1956). Background 
information for algorithms that use continued fractions may be found in H. S. Wall’s Analytic Theory of 
Continued Fractions (D. VanNostrand Co., Inc., Princeton, N. J., 1948). 


IHCCLABS (CDABS) and IHCCSABS (CABS) Subprograms 
1. Write |x + iy| = a + ib. 
2. Ifx = y = 0, thena = 0 and b = 0. 
3. Let vo, = max (|x|, [y|), and 
02 = min (|x|, |y]). 


= fig Ne. 
Then, a = oa] 1+ (=) , and b = 0. 
1 
The algorithms for both complex absolute value subprograms are identical. 


Each subprogram uses the appropriate real square root subprogram (IHCLSQRT 
Or IHCSSQRT ). 





IHCCLEXP (CDEXP) and IHCCSEXP (CEXP) Subprograms 


Algorithm 


The value of e?+” is computed as e” * cos(y) + i* e” * sin(y). The algorithms 
for both complex exponential subprograms are identical. Each subprogram uses 
the appropriate real exponential subprogram (IHCLEXP or IHCSEXP ) and the appro- 
priate real sine/ cosine subprogram (THCLSCN or IHCSSCN ). 


Effect of an Argument Error 


The effect of an argument error depends upon the accuracy of the individual 
parts of the argument. If e7+” = R « e'#, then H = y and e(R) ~ A (x). 


IHCCLLOG (CDLOG) and IHCCSLOG (CLOG) Subprograms 


Algorithm 
1. Write log, (x + iy) =a + ib. 
2. Then, a = log, |x + iy|and b = the principle value of arctan 4 : 


The algorithms for both complex natural logarithm subprograms are identical. 
Each subprogram uses the appropriate complex absolute value subprogram 
(IHCCLABS or IHCCSABS), the appropriate real natural logarithm subprogram 
(IHCLLOG or IHCsLOG), and the appropriate arctangent subprogram (IHCLATN2 
Or IHCSATN2). 


Effect of an Argument Error 


The effect of an argument error depends upon the accuracy of the individual 
parts of the argument. If x + iy = r¢ e* and log, (x + iy) = a + ib, thenh = b 
and E(a) = 8(r). 


IHCCLSQT (CDSQRT) and IHCCSSQT (CSQRT) Subprograms 
Algorithm | 

1. Write Vx + iy = a + ib. 

2. Ifx = y = 0, thena = Oandb = 0. 


3. Ife 0,thena = | 2 - Et iy| 
andb = ¥_, 
2a 
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4, Ifx < 0, thenag = 
2a 


db= e a. 
an (sign y) N 


The algorithms for both complex square root subprograms are identical. Each 
subprogram uses the appropriate real square root subprogram (IHCLSQRT or 
THCSSQRT ). 


Effect of an Argument Error 
The effect of an argument error depends upon the accuracy of the individual 
parts of the argument. If x + iy = r¢ e” and Vx + iy = R° e, then 


nee = (1), and ¢(H) anche): 


IHCCLSCN Subprogram (CDSIN and CDCOS) 
Algorithm 


l. If the sine is desired, then 
_sin(x + iy) = sin(x) * cosh(y) + i* cos(x) * sinh(y). 
If the cosine is desired, then 
cos(x + iy) = cos(x) * cosh(y) + i* sin(x) * sinh(y). 


2. Ifx < 0, then sinh( —x) = —sinh(<x). 


3. Ifx > 0.3465736, then sinh(x) = 





2 


4. If0 S x S 0.3465736, then compute sinh (x) by use of the polynomial: 
— = Ay + ax? + doxt +... + asx?. 
The coefficients are obtained by expanding the polynomial with respect to the 
Chebyshev polynomials over the range 0 = x? = 0.120113. The relative error 
of this approximation is less than 2—21-8, 
5. The value of cosh(x) is computed as cosh(x) = sinh|x| + a ; 
e 
This computation uses the real exponential subprogram (1HCLExP) and the 
real sine/cosine subprogram (IHCLSCN). 


Effect of an Argument Error 

To understand the effect of an argument error upon the accuracy of the answer, 
the programmer must understand the effect of an argument error in the mHCLSCN, 
IHCLEXP, and IHCLSNH subprograms. 


IHCCSSCN Subprogram (CSIN and CCOS) 


Algorithm 
1. If the sine is desired, then 
sin (x + iy) = sin(x) * cosh(y) + i* cos(x) * sinh(y). 
If the cosine is desired, then 
cos(x + iy) = cos(x) * cosh(y) + i* sin(x) * sinh(y). 
2. Ifx < 0, then sinh( —x) = —sinh(x). 


3. Ifx > 0.3465736, then sinh(x) = ——— .. 


4. If0 =x <= 0.3465736, then compute sinh (x) by use of the polynomial: 
sinh(x) 
x 


= dy + ayx? + aox'. 


The coefficients are obtained by expanding the polynomial with respect to the 
Chebyshev polynomials over the range 0 = x? = 0.120113. The relative error 
of this approximation is less than 2-764, 


5. The value of cosh(x) is computed as cosh(x) = sinh|x| + = , 
e OH 


This computation uses the real exponential subprogram (1HcsExp) and the 
rea] sine/cosine subprogram (IHCSSCN ). 


Effect of an Argument Error 


To understand the effect of an argument error upon the accuracy of the answer, 
the programmer must understand the effect of an argument in the rHcsscn, 
THCSEXP, and IHCSsCNH subprograms. 


IHCLASCN Subprogram (DARSIN and DARCOS) 


Algorithm 
1 f0O<*s + , then compute arccos(x) as: 
arccos(x) = igo arcsin(x). 


fos*s es then compute arcsin(x) by a polynomial of the form: 


— 2 > 
arcsin(x) = x + €4x® + cox5 +... Cy2x?5, 
i : : arcsin(x) 
The coefficients are obtained by expanding the function f(z) = =p 
1 


% = x*, with respect to the Chebyshev polynomials over the range, 0 = x ST: 


The relative error of this approximation is less than 2-557, 
25 Te >< x <= 1, then compute arcsin (x) as: 


arcsin(x) = oo arccos(x). 
If + <x <1, then compute arccos (x) as: 


arccos(x) = 2° arcsin (J : 7 *) : 


This case is now reduced to the first case because within these limits, 
1-«x < 1 : 
2. =. 9 








0S 


This computation uses the real square root subprogram (IHCLSQRT). 
3. If -1Sx < 0, then arcsin(x) = —arcsin|x| 
and arccos(x) = 2 — arccos|x|, 


This reduces these cases to one of the two positive cases. 


Effect of an Argument Error 
tA 
V1l— x? 
a small A causes a substantial error in the answer. For the arcsine, « ~ 8 if the 
value of x is small. 


Rex 





. For small values of x, E ~ A. Toward the limits (+1) of the range 


Appendix A. Algorithms 21 


IHCLATAN Subprogram (DATAN) 


Algorithm 
1, Reduce the computation of arctan (x) to the case 0 = x = 1 by using 
arctan( —x) = —arctan(x) or 
1 7 
arctan a a arctan |x| 


2. If necessary, reduce the computation further to the case |x| = tan 15° by using 


arctan (x) = 30° + arctan eA seat =i 
x+ V3 
V3e°x—1 
x+ V3 
tan 15° <x1. The value of (\/3*x— 1) is computed as 
1 


(V3 —-1)x- > = + x to avoid the loss of significant digits. 


3. For |x| S tan 15°, use a continued fraction of the form: 


The value of <= tan 15° if the value of x is within the range, 








arctan (x) Agee a,x? dy az a4 
x re (by + x?) + (be + x?) + (bs + x?) + (bg + x?) + °°" 


The relative error of this approximation is less than 2—5*°. The coefficients of 
this formula were derived by transforming the continued fraction: 





wk. 3°4 16 * 25 
arctan (x) _ 1+ 3 25° 7 7° 81°11 
% i. Soy y-2) 23 + x-2)} — DO te tid = 
5 5°9 9°13 
4*3°9 
5° 11+ 169 





can 
(Gaea7 ** ) 


9 
‘ ee eee mee 
where w has an approximate value of Bell 13°17 (— x —2 + 40) but the true 
64 + 27 
5+ 289 + 19 


value of w is 179 wees 
poe ieie netrea —2 
(a-7: ne )+ 


Effect of an Argument Error 
A 

l+x° 

effect of « upon § diminishes. 


For small values of x, « ~ 8, and as the value of x increases, the 


IHCLATN2 Subprogram (DATAN and DATAN2) 
Algorithm 


v 


7° 


x1 


1. For arctan(%;, x2), if either x. = 0 or 











> 258 the answer = (sign x,) ° 
Otherwise, if x. > 0, the answer = arctan ( =) , and 
2 


if x2 < 0, the answer = arctan ( = ) + (sign x,) * x. 
2 


The rest of the computation is identical for either one or two arguments. 


22 


. Reduce the computation of arctan(x) to the case 0S x S 1, by using 


arctan( —x) = —arctan(x), or 


arctan Le) ee arctan|x| 
rs 3 


. If necessary, reduce the computation further to the case |x} < tan 15° by using 


arctan(x) = 30° + arctan ae Soe 
x+ V3 
V3°x-1 oa che edie 
The value of |-———————} S tan 15° if the value x is within the range, 
x+ V3 








tan 15° < x <1. The value of (1/3 * x — 1) is computed as (V3 — 1) x—1 
to avoid the loss of significant digits. 


. For |x| <S tan 15°, use a continued fraction of the form: 


arctan (x) ee ax? dy as a4 
« =. (by +x?) + (be +x?) + (bs + x7) + (bg + x?) + 


The relative error of this approximation is less than 2—57-°, The coefficients of 
this formula were derived by transforming the continued fraction: 
= = 3°4 16° 25 
arctan(x) _ 95°7 7°8l°11 


x Ss ae = 39 -2\_ 
“Gea Gry Ge 
4*3°9 
5° 11¢ 169 


179 a 
(Ga Wt * ) Ww 











: 2 sats 
where w has an approximate value of 5-1 13°17 (—x-2 + 40) but the 
64 * 27 
5 * 289 + 19 


true value of w is 


WO ya ot 
(seat *)+ 


Effect of an Argument Error 


A . 
E ~ ——.. For small values of x, « ~ 8, and as the value of x increases, the effect 


] + x? 


of « upon § diminishes. 


IHCLERF Subprogram (DERF and DERFC) 
Algorithm 


1. 


If 0 =x < 1, then compute the error function by the following approximation: 
erf(x) Sx (dy + ayx? + agxt +... + ayyx?*), 

The coefficients were obtained by expanding the function f(z) = one), 

z = x?, with respect to the Chebyshev polynomials over the range, 0 =x < 1. 

The relative error of this approximation is less than 1.07 » 2~57, The value of 

the complemented error function is computed as erfc(x) = 1 — erf(x) and is 


greater than <= ; 


. If l= x < 2.0400009, then compute the complemented error function by the 


following approximation: 
erfc (x) > bo + biz + bez? = ae ees 5 by 5218 
where z = x — Ty and Ty = 1.999999,¢. The coefficients were obtained by ex- 
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panding the function f(z) = erfc(z + T)) with respect to the Chebyshev poly- 
nomials over the range —1 = x = 0.04. The absolute error of this approxima- 
tion is less than 1.5 * 2—®!, The limits of this range and the base value for To 
were used to minimize the hexadecimal truncation error. The value of the com- 


1 
plemented error function within this range is greater than = . The value of 


256 
the error function is computed as erf(x) = 1 — erfe(x). 


3. If 2.0400009 = x = 13.306, then compute the complemented error function by 
the following approximation: 


Co + Cyx7? + Cox 4 He... HK Capx?9) e-” 
A ee ec ae ce ce : = = au 


The coefficients were obtained by expanding the function 
f(z) = erfc (x) * x * e”, z = x~?, with respect to the Chebyshev polynomials 
over the range 2.04-2 > z = 13.3062. The relative error of this approximation 
ranges from 2-58 at 2.04 to 2-5! at 13.306. This computation uses the real ex- 
ponential subprogram (IHCLEXpP ). 
If x <= 6.092, then the error function is computed as erf(x) = 1 — erfc(x). 
If x > 6.092, then the error function is = I. 
4, If 13.306 < x, then the error function is = 1, and the complemented error 
function is = 0. 
5. If x < 0, then reduce to a case involving a positive argument by the use of 
the following formulas: 
erf(—x) = —erf(x) and erfe(—x) = 2 — erfc(x). 
Effect of an Argument Error 
E ~ e—” A. For the error function, as the magnitude of the argument exceeds 
1, the effect of an argument error upon the final accuracy diminishes rapidly. For 


small values of x, « ~ 8. For the complemented error function, if the value of 


x is greater than 1, erfe(x) ~ a. Therefore, « ~ 2x° * 8. If the value of x 
x 


is negative or less than 1, thene ~ e-”° A, 


IHCLEXP Subprogram (DEXP) 

Algorithm 

1. If x << — 180, 2183, then 0 is given as the answer. 
2. Divide x by log,2 and write 


x Cc 
= = (4a — b -—— -d 
m log,2 ( 16 ) 

where a, b, and ¢ are integers, OS b<3, 0c 15, and d is within the range 


OSd< > . Then, e? = 2% = 164° 2-2 » 2-0/16 « 9—4, 





3. Compute 2—¢ by using the Chebyshev interpolation of degree 6 over the range, 


Osde< — . The maximum relative error of this computation is 2-5". 


4, Ifc > 0, then multiply 2-4 by 2-°/1%, (The 15 values of 2-°/"6 forl Sc 15 
are included in the subprogram. ) 

5. If b > 0, then halve the result b times. 

6. Finally, add the hexadecimal exponent a to the characteristic of the answer. 


Effect of an Argument Error 


E ~ A. If the magnitude of x is large, even the round-off error of the argument 
causes a substantial relative error in the answer because A = 8° x. 


IHCLGAMA Subprogram (DGAMMA and DLGAMA) 


Algorithm 

1, If 0 < x S 2-52, then compute log-gamma as log.I'(x) = —log.(x). This 
computation uses the real logarithm subprogram (IHCLLOG). 

2. If 2-252 < x < 8, then compute log-gamma by taking the natural logarithm 
of the value obtained for gamma. The computation of gamma depends upon 
the range into which the argument falls. 


3. If 2-952 < x < 1, then use I(x) = Mat’) to reduce to the next case. 


4. If 1 = x S 2, then compute gamma by the following approximation: 
T(x) = do + 4% + dee +... HF deez” 

where z = x — 1.5. The coefficients were obtained by expanding the function 
f(z) = I(x) with respect to the Chebyshev polynomial for |z| = 0.5. The 
absolute error of this approximation is less than 1.5 + 2-58, 

5. If2 <x < 8, then use T(x) = (x — 1) I(x — 1) to reduce to the preced- 
ing case. 

6. If 8 = x, then compute log-gamma by the use of Stirling’s formula: 


log.l (x) = x(loge(x) — 1) — s log,(x) + > loge(Qn) + G(x). 


The modifier term G(x) is computed as 
G(x) = byx-! + box? + bex 5 + bax-t + bdsx-. 

The coefficients were obtained by expanding the function f(z) = Cts) ; 
z = x, with respect to the Chebyshev polynomials over the range 
0 < z < 8-*. The absolute error of the approximation for G(x) is less than 
x * 2-58 Because, in this range, x < log.(x), the contribution of this error 
to the relative error of the value for log-gamma is less than 2—5*. This com- 
putation uses the real logarithm subprogram (1HCLLOG). 

For. gamma, compute [(x) = e”, where y is the value obtained for log- 


gamma. This computation uses the real exponential subprogram (IHCLEXxP). 





Effect of an Argument Error 
e ~ w(x) * A for gamma, and E ~ y(x) * A for log-gamma, where y is the 
digamma function. 


If 5 <x < 38, then —2 < y(x) < 1. Therefore, E ~ A for log-gamma. 


However, because x = 1 and x = 2 are zeros of the log-gamma function, even 
a small § can cause a substantial « in this range. 

If the value of x is large, then y(x) ~ log.(x). Therefore, for gamma, 
e ~ 8° x ° log.(x). In this case, even the round-off error of the argument con- 
tributes greatly to the relative error of the answer. For log-gamma with large 
values of x, « ~ 5. 


IHCLLOG Subprogram (DLOG and DLOG10) 


Algorithm 
1. Write x = 16? + 2-¢+ m, where p is the exponent, q is an integer, 0 q $3, 


and m is within the range, = =m<l 
2. Define two constants, a and b (where a = base point and 2~» = a) as follows: 
] i} 1 
— <= —_—— = — — 
If 7 am< wan 5 and b = 1. 
If agen <1, thena = landb = 0. 
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6. 


. Write z = 
. Now, x = 247-0 - 


. Finally, log. ( 


m—-a 2 ta ae 
mig? phen. m =a Toe 








and |z| < 0.1716. 


+ z 
—Zz 





Sndioped (an 9th seo iow: (724). 


; a is computed by using the Chebyshev interpolation of 
—Z 

degree 7 in z? over the range, 0 z? = 0. 02944. The maximum relative error 
of this approximation is 2—5°-6, 


If the common logarithm is desired, then logiox = logioe * log.x. 





Effect of an Argument Error 
E ~ 8. Therefore, if the value of the argument is close to 1, the relative error can 
be very large because the value of the function is very small. 


IHCLSCN Subprogram (DSIN and DCOS) 
Algorithm 


1. 


3. 


. Finally, either sin 


Divide |x| by oe and separate the quotient (z) into its integer part (q) and 


ee q + 1, where q is an integer and r 
TT 





its fraction part (r). Then, z = |x 





is within the range, Or < 1. 


. If the cosine is desired, add 2 to q. If the sine is desired and if x is negative, 


add 4 to q. This adjustment of q reduces the general case to the computation of 
sin (x) for x = 0, because 


cos (+x) = sin ( ++) , and 


sin (—x) = sin (|x| + 7). 
Let qo = q mod 8. 


Then, for go = 0, sin (x) = sin (+ r) 

go = 1, sin (x) = cos (Fa _ r)) 

do = 2, sin (x) = cos ( +r) 

qo = 3,sin (x) = sin (4a = ")) 
eo) 
do = 5, sin (x) = —cos (Fa - r)) 
do = 6, sin (x) = —cos (+7) 


do = 7, sin (x) = —sin (Fa = r)) 


These formulas reduce each case to the computation of either sin (= r) or 








do = 4, sin (x) = —sin ( 





cos * 11 }; where 7, is either r or (1 --r), and is within the range, O71 =1. 


(< 


( i * 71 } or cos i 1, ) is computed, using the Chebyshev 


interpolation of degree 6 in 1,? for the sine, and of degree 7 in r,? for the cosine. 
The maximum relative error of the sine polynomial is 2—5* and that of the 
cosine polynomial is 2— 6-8, 





Effect of an Argument Error 

E ~ A. As the value of the argument increases, A increases. Because the function 
value diminishes periodically, no consistent relative error control can be main- 
T 


2° 


Tv 


<xy< 
5) pay i 


tained outside of the principal range, — 


IHCLSCNH Subprogram (DSINH and DCOSH) 


Algorithm 


1. If |x| < 0.3465736, then compute sinh(x) as: 
sinh(x) = x + cyx® + cox® + Cgx™ + c4x® + cx}. 
The coefficients are obtained by expanding the function f(z) = sini) ; 
x 
z = x*, with respect to the Chebyshev polynomials over the range, 


0 =z < 0.12011326, The relative error of this approximation is less than 2—°1-°, 
2. If either |x| 2 0.3465736 or the cosh(x) is desired, obtain w = e!*!. Then, 
w+w7t =i 


cosh(x) = — 7 and sinh(x) = (sign x) = 5 


nential subprogram (IHCLEXP) is used to compute the value of w. 


1 
. The real expo- 


Effect of an Argument Error 


For the hyperbolic sine, E ~ A * cosh(x) and e ~ A * coth(x). 
For the hyperbolic cosine, E ~ A sinh(x) and «e ~ A ®* tanh(x). 

Specifically, for the cosine, E ~ A over the entire range; for the sine, « ~ 8 
for the small values of x. 


IHCLSQRT Subprogram (DSQRT) 


Algorithm 
1. If x = 0, then the answer is 0. 


2. Write x = 167?—2* m, where 2p — q is the exponent and q equals either 0 or 1; 
m is the mantissa and is within the range, = =Sm=l. 
3. Then, Vx = 16? * 2-24+ \/m. 


4, For the first approximation of \/x , compute the following: 


pa 
yo = 2% 16r+ (5-4 9 ‘m). 





The maximum relative error of this approximation is 7° 


5. Apply the Newton-Raphson iteration 


1 
rs (08) 


four times to yo (the first two times in the short form and the last two times in 
the long form). The final step is performed as 


_ lfx 
Ys = Ys t 7a ere vs) 
to minimize the computational truncation error. The maximum relative error 
of the final result is theoretically 2—%-”°, 


Effect of an Argument Error 


1 


e~ —9 
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IHCLTANH Subprogram (DTANH) 
Algorithm 


1. If |x| < 0.54931, then use the following fractional approximation: 


tanh (x) Pa ee 4x2 + doxt + agx® + x8 
x bo + byx? + boxt + bgx® + x8 


where: 
a, = 676440.765 bo = 2029322.295 
dz = 45092.124 b, = 947005.29 
d3 = 594.459 be = 52028.55 
bz = 630.476 


The maximum relative error of this approximation is 2—°-5, The formula was 
obtained by transforming the continued fraction 
tanh (x) _ yy 2 x 
83+ 5+°°' 15+w 
where w has an approximate value of 0.017, but the true value of w is 
x x 
17+ 19+ °°" 
2. If 0.54931 = x < 20.101, then use the identity tanh (x) = 1 — Ta 


computation uses the double precision exponential subprogram (IHCLEXP). 
3. Ifx = 20.101, then tanh (x) = 1. 
4, Ifx <= — 0.54931, then use the identity tanh (x) = — tanh ( —x). 





. This 


Effect of an Argument Error 


E ~ (1 — tanh? x) A, ande ~ . For small values of x , « ~ 8. As the 


2A 
sinh (2x) 
value of x increases, the effect of § upon ¢ diminishes. 


IHCLTNCT Subprogram (DTAN and DCOTAN) 


Algorithm 


1. Divide |x| byt and separate the result into the integer part (q) and the 
fraction part (1). Then, |x| = rn (q + 1). 


2. Obtain the reduced argument (w) as follows: 
if q is even, then w = 1. 
if q is odd, then w = 1 — 1, 
The range of the reduced argument is 0 = w = 1. 


3. Let go = q mod 4. 
Then, for go = 0, tan |x| = tan (+ °w ) and cot |x| = cot (4 ° w) 
do = I, tan |x| = cot ( a “w ) and cot |x| = tan (4 -w) 
go = 2, tan |x| = —cot ( : ° w) and cot |x| = tan ( 7 -w) 


qo = 3, tan |x| = tan ( Zi -w) and cot |x| = cot ( rk -w) 




















4. The values of tan oe w) and cot ( i ° w) are computed as the ratio 


of two polynomials. 








tan Fw) = ory ace ‘w) =o 
we P(w*) 


where P(w?) is of degree 3 and Q(w?) is of degree 4 in w?. The coefficients 
of P and Q are obtained by economizing the continued fraction 


tan(z) _ on zm 
z a. 52: 7 
in the following way. 
2 2 2 2 2 9 2 
Write: @2(2) 27 2 #@ #@ # _ % z z 


and determine the values for d,, da, and ds so that the right-hand expression 
gives the exact answers for z? = 0.395, 0.542, and 0.607. Then the maximum 


relative error of this formula over the range 0 = z S rs is 3.4 * 10-19, 
Change the variable from z to w = 7 * z and rewrite the formula to obtain 
P(w?) and QO(w). 

5. If x < 0, then tan(x) = —tan 





, and cot(x) = —cot |x|. 


Effect of an Argument Error 


2 ; ya 
~ =z» and « ~ ——— for tan(x). Therefore, near the singularities 
cos? (x) sin{ 2x) 
1 eer 
ofx = k++ > where k is an integer, no error control can be maintained. 


This is also true for cotan(x) for values of x near kz, where k is an integer. 


IHCSASCN Subprogram (ARSIN and ARCOS) 


Algorithm 


Lif0sxcs + then compute arccos(x) as 
arccos(x) = oe arcsin(x). 
[foOs*cs > then compute arcsin(x) by a polynomial of the form: 
arcsin(x) sx + 04x? + Cox® + cgx™ + c4x® + osx. 
; : arcsin(x ) 
The coefficients are obtained by expanding the function f(z) = oe rs 
z = x, with respect to the Chebyshev polynomials over the range 0 <= z <-. 


The relative error of this approximation is less than 2—27-5, 
2. Tf + < x <= 1, then compute arcsin(x) as: 
arcsin(x) = 5 — arccos(x). 
If + < x 1, then compute arccos(x) as: 


arccos(x) = 2° aresin (J ). 


This case is now reduced to the first case because within these limits, 





This computation uses the real square root subprogram (IHCSSQRT). 


3. If -1 Sx < 0, then aresin(x) = — arcsin |x|, and arccos(x) = 7 — arccos |x|. 
This reduces these cases to one of the two positive cases. 
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Effect of an Argument Error 


Ries ee oe small values of x, E ~ A. Toward the limits (+1) of the 


V1l—+# 


range, a small A causes a substantial error in the answer. 


IHCSATAN Subprogram (ATAN) 


Algorithm 
1, Reduce the computation of arctan (x) to the case 0 = x = 1, by using 


arctan (—x) = —arctan (x), or 
arctan ip eae arctan |x|. 
Ix] J 2 
2. If necessary, reduce the computation further to the case |x| <= tan 15° by using 


arctan (x) = 30° + arctan (oS — *). 


_ + V3 

io ta 1 ey OB eid 

—————-| <= tan 15° if the value of x is within the range, 
x+ V3 

tan 15° < x <1. The value of (1/3 * x — 1) is computed as 

(V3 —1) x — 1 + x to avoid the loss of significant digits. 


3. For |x| S tan 15°, use the approximation formula: 


arctan (x) __ - 0.55913709 
a 0.60310579 — 0.05160454x? + + 3 14087812 +. 14087812 


This formula has a relative error less than 2-271 and can be obtained by 
transforming the continued fraction 


The value of 








x2 
arctan (x) _ l en 5 
x 3+ ( 5 *) 
a —w 
v4 
where w has an approximate value of ( — Be —2 4 ste 0-4, but the true 
4°5 
7*°7*9 


value of w is —————————_ 


ee aicrt em —2 
(3 FI ee 3 )+ 
The original continued fraction can be obtained by transforming the Taylor 


series into continued fraction form. 


Effect of an Argument Error 


A ‘ 
Ew Tae - For small values of x, « ~ 8; as the value of x increases, the effect 


of § upon e diminishes. 


IHCSATN2 Subprogram (ATAN and ATAN2) 
Algorithm 
> 274, the answer = (sign x1) * 








1. For arctan(x;, x2), if either x. = 0 or |— : oe 


2 


: . x: 
Otherwise, if x. > 0, the answer = arctan =) , and 
2 


if x2 < 0, the answer = arctan (+) + (sign x,) 7. 
2 


The rest of the computation is identical for either one or two arguments. 


2. Reduce the computation of arctan(x) to the case 0 S x = 1, by using 
arctan( —x) = — arctan(x) or 


l1\_@- 
arctan (a)= 5} arctan |x|. 


3. If necessary, reduce the computation further to the case |x| S tan 15° by using 


arctan(x) = 30° + arctan Motes 
s+ V3 
V3ex-—1 ; : eae 
The value of |>————| SS tan 15° if the value of x is within the range, 
xt V3 








tan 15° < x S 1. The value of (3 ° x — 1) is computed as 
(V3 — 1) x — 1+ x to avoid the loss of significant digits. 


4. For |x| tan 15°, use the approximation formula: 
arctan(x) | a2 , 0.55913709 
ao 0.060310579 — 0.0516045x? + x + 14087813 
This formula has a relative error less than 2~27:! and can be obtained by trans- 
forming the continued fraction 


x2 
arctan(x) 8 
x = 3+ ( 5 a) 
a ee = Ap 
7 
where w has an approximate value of ( - Pe + =) 10-4 but the true 
4°5 
7°7°9 
value of w is 773 — 2 
paeeer nies —2 
(Fat? )+ | 


The original continued fraction can be obtained by transforming the Taylor 
series into continued fraction form. 


Effect of an Argument Error 


A : 
Ew Tee: For small values of x, « ~ 5; as the value of x increases, the effect 


of § upon e diminishes. 


IHCSERF Subprogram (ERF and ERFC) 
Algorithm 


1. If 0 =x S 1.317, then compute the error function by the following approxi- 
mation: 
erf(x) 2 x(do + ax? + doxt +... + agx!?), 


The coefficients were obtained by expanding the function f(z) = eae oe 


erf (x) 
x 
with respect to the Chebyshev polynomials over the range 0 = x = 1.317. The 
relative error of this approximation is less than 2-24, The value of the comple- 


mented error function is computed as erfc(x) = 1 — erf(x) and is greater 


1 
than T6" 


2. If 1.317 < x & 2.0400009, then compute the complemented error function by 
the following approximation: 
erfe(x) = bo + byz + baz? +... + byt 
where z = x — To and Ty = 2.0400009. The coefficients were obtained by 
expanding the function f(z) = erfc(x + To) with respect to the Chebyshev 
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polynomials over the range (1.317 — Ty) < z S 0. The absolute error of this 
approximation is less than 1.3 * 2—9°, The value of the complemented error 
function within the range 1.317 < x <= Tp is greater than = . The value of the 
error function is computed as erf(x) = 1 — erfe(x). 


3. If To < x = 13.306, then compute the complemented error function by the fol- 

lowing approximation: 
(Co + cyx—2 + cox—4* +... + Cgx—1?) e-” 

x 

The coefficients were obtained by expanding the function 
f(z) = erfe (x) + x* e”, z = x—2, with respect to the Chebyshev polynomials 
over the range T)-? > z = 13.306—?. The relative error of this approximation 
is less than 1.2 * 2-28, This computation uses the real exponential subprogram 
(THCSEXP ). 
If x = 3.9192, then the error function is computed as erf(x) = 1 — erfc(x). 
If x > 3.9192, then the error function is = 1. 


4. If 13.306 < x, then the error function is =< 1, and the complemented error func- 
tion is = 0, 


erfc(x) = 


5. If x < 0, then reduce to a case involving a positive argument by the use of the 
following formulas: 
erf(—x) = —erf(x) and 
erfe( —x) = 2 — erfc(x). 


Effect of an Argument Error 


E ~ e—” A. For the error function, as the magnitude of the argument exceeds I, 
the effect of an argument error upon the final accuracy diminishes rapidly. For 
small values of x, « ~ 8. For the complemented error function, if the value of x is 


greater than 1, erfc(x) ~ <— . Therefore, « ~ 2 x? + 8, If the value of x is negative 


or less than 1, then e ~ e-”° A. 


IHCSEXP Subprogram (EXP) 
Algorithm 
1. Ifx < —180.218, then 0 is given as the answer. 
2. If |x| < 2—?8, then 1 is given as the answer. 
3. Otherwise, divide x by log.2 and write 
aad a ia (4a — b - d) 
where a and b are integers, OS bS3and0Sd<1. 
Then, e” = 24 = (164+ 2-%+2-¢), 
4. Compute 2—4 by the following fractional approximation: 
2-4 = cats 


0.034657359 d? + d + 9,9545058 — _ 0L1.97227 


d? + 87,417497 


This formula can be obtained by transforming the Gaussian-type continued 
fraction 


1+ 2— 3+ 2—5+ 2—7+ 2 
The maximum relative error of this approximation is 2—°. 
5. Multiply 2-4 by 2-?. 


6. Finally, add the hexadecimal exponent a to the characteristic of the answer. 


Effect of an Argument Error 
« ~ A. If the magnitude of x is large, even the round-off error of the argument 
causes a substantial relative error in the answer because A = 8 ° x. 


IHCSGAMA Subprogram (GAMMA and ALGAMA) 

Algorithm 

1. If0 < x S 2-2, then compute log-gamma as log, T (x) = —log, (x). This 
computation uses the real logarithm subprogram (IHCSLOG ). 

2. If 2-252 < x < 8, then compute log-gamma by taking the natural logarithm of 
the value obtained for gamma. The computation of gamma depends upon the 
range into which the argument falls. 


(x +1) 
x 


3. If2-282 < x < 1, then useI (x) = x to reduce to the next case. 


4. If 1 Sx S 2, then compute gamma by the following approximation: 
T(x) = do + a,% + 92” +... + doz? 
where z = x — 1.5. The coefficients were obtained by expanding the function 
f(z) = I’ (x) with respect to the Chebyshev polynomials for |z| <= 0.5. The 
absolute error of this approximation is less than 1.5 * 2-5, 
5. If2 <x < 8, then use (x) = (x — 1) T (x — 1) to reduce step by step to 
the preceding case. 
6. If 8 < x, then compute log-gamma by the use of Stirling’s formula: 
loge F (x) 2 x(loge (x) — 1) — slog, (x) +5 loge (2r) + G(x). 
The modifier term G(x) is computed as 
G(x) = byx-1 + box-?. | 
The absolute error of the approximation for G(x) is 1.4 * 2-28. This computa- 
tion uses the real logarithm subprogram (1HCSLOG). 
For gamma, compute [ (x) = e”, where y is the value obtained for log- 
gamma. This computation uses the real exponential subprogram (IHCSEXP ). 


Effect of an Argument Error 
e ~ y (x) * A for gamma, and E ~ y (x) « A for log-gamma, where y is the 
digamma function. 


If = <x <3, then —2 < y (x) < I. Therefore, E ~ A for log-gamma. How- 


ever, because x = 1 and x = 2 are zeros of the log-gamma function, even a small 
8 can cause a substantial ¢ in this range. 

If the value of x is large, then y (x) ~ log. (x). Therefore, for gamma, 
e ~ 8x ° log, (x). In this case, even the round-off error of the argument con- 
tributes greatly to the relative error of the answer. For log-gamma with large 
values of x, « ~ 6. 


IHCSLOG Subprogram (ALOG and ALOG10) 

Algorithm 

1. Write x = 16? + m, where p is an integer and m is within the range,+5 sm<l 
2. Define two constants, a and b, where a = base point and 2—” = a, as follows: 


It =m< =, then a Soha = 4, 


16 
1 1 =, 2h _ 
tf Sm < > thena = Z and b = 2. 
If-5- Sm <1,thena =1andb = 0, 
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m—a l+2z 1 
i = = e es <= 
3. Write z = . Then, m = a ee ,and |z| <= z 
4. Now, x = 2479-? « ; : , and log.x = (4p — b) log.2 + logs ( ' = ). 


: = *) is evaluated using the Chebyshev interpolation of degree 





5. Finally, log. ( 


: 1 ; ; ; 
4 in z? over the range, OS 2?7= >: The maximum relative error of this ap- 
proximation is 2—27-8, 


6. If the common logarithm is desired, then logiox = logice * log.x. 


Effect of an Argument Error 


E ~ 68. Specifically, if 8 is the round-off error of the argument, e.g., 8 ~ 6 * 10-8, 
then E ~ 6° 10-8, Therefore, if the argument is close to 1, the relative error can 
be very large because the value of the function is very small. 


IHCSSCN Subprogram (SIN and COS) 
Algorithm 


l. Define z = ae 
T 


(7). Then z = q +1, and |x| = (7:7) +(4°). 


2. If the cosine is desired, add 2 to q. If the sine is desired and if x is negative, 
add 4 to q. This adjustment of q reduces the general case to the computation 
of sin (x) for x = 0 because 


cos (+x) = sin (++ *) , and 





x| and separate z into its integer part (q) and its fraction part 





2 
sin(—x) = sin(r +x). 
3. Let qo = q mod 8. 
Then, for qo = 0, sin (x) = sin ( = r) 





do = 1, sin (x) = cos 





qo = 3, sin (x) = sin 


( 
do = 2, sin (x) = cos ri 
(+ 





4 
T 
do = 4, sin (x) = sin ( ° ) 


qo = 5, sin (x) = ~cos ( 
do = 6, sin (x) = -008 (7 -r) 


do = 7, sin (x) = —sin (t=) 








These formulas reduce each case to the computation of either sin ( 7 . ri) or 


cos (} 1 ) where 1; is either 7 or (1 — 1) and is within the range, 


0Ofn=1. 


4, Finally, the computation for either the sine or the cosine is performed, using 
the Chebyshev interpolation of degree 3 in 7,7. The maximum relative error 
of the sine polynomial is 2—28! and that of the cosine polynomial is 2—*4-°. 


Effect of an Argument Error 


E ~ A. As the value of x increases, A increases. Because the function value dimin- 

ishes periodically, no consistent relative error control can be maintained outside 
. e TT Tv 

the principal range, — 7= xs+ 7° 


IHCSSCNH Subprogram (SINH and COSH) 


Algorithm 
1. If |x| < 0.3465736, then compute sinh(x) as: 
sinh(x) = x + 0.16666505x? + 0.00836915x5, 


The coefficients were obtained by expanding the function f(z) = ts) ’ 


= x, with respect to the Chebyshev polynomials over the range 
< z < 0.12011326. The relative error of this approximation is less than 
— 26.5 


NON 


2. If either |x| 2 0.3465736 or the cosh(x) is desired, obtain w = e'!, Then, 
=a _ 1 
oe, and sinh(x) = (sign x). — The real exponen- 


tial subprogram (1HCSEXxP) is used to compute the value of w. 


cosh(x) = 


Effect of an Argument Error 


For the hyperbolic sine, E ~ A* cosh(x) and e ~ A» coth (x). 
For the hyperbolic cosine, E ~ A* sinh (x) ande ~ 8° tanh(x). 

Specifically, for the cosine, « ~ A over the entire range; for the sine, e ~ 8 for 
small values of x. | 


IHCSSQRT Subprogram (SQRT) 
Algorithm 
1. Ifx = 0, then the answer is 0. 


2. Write x = 16? * m, where p is an integer and m is within the range, 


1 
3. Then, \/x = 16? + \/m, where p is the exponent of the answer and m is the 
mantissa of the answer. 


4. For the first approximation of \/m, take hyperbolic approximations of the form 





at ; 2 ; where the values of a, b, and c depend upon the value of m as follows: 
a. ts <m <1,then a = 180713 
= —1,57727 
c = 0.954182 


_ These values minimize the maximum relative error (9) over the range, 
while making an exact fit at m = 1. The exact fit at m = 1 minimizes the 
computational loss of the last hexadecimal digit for the values of m slightly 
less than 1. The relative error of this approximation is less than 2—5-44, 
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b. If = an <ag ‘then a = 0.428795 
= —0,0214398 
c = 0.0548470 


These values minimize m1/8 * «9 over this range of m where «9 denotes the 
relative error of this approximation. ¢o is less than 2—&5 « m-1/8, 
5. Multiply the result by 16? to obtain the first approximation (yo) of the answer. 
6. To obtain the final answer, the Newton-Raphson iteration 


= x 


must be applied twice to yo. For 16 <= m <1, the final relative error is 


theoretically less than 2—74:7; for 7 =m <@ the final absolute error is 


theoretically less than 2—?° + 16?. 


Effect of an Argument Error 


e~Z 8. 


2 


IHCSTANH Subprogram (TANH) 


Algorithm 
1. If |x| = 2-™, then tanh (x) = x. 
2. If 2-12 < |x| < 0.54931, use the following fractional approximation: 
tanh (x) ae x? + 35.1535 
: x? + 45,1842 +P? 
This approximation has a relative error less thari 2—2*. The formula can be 
obtained by transforming the continued fraction 





tanh (x) aa Me a 
x 3+ 5+ 7+w 
where w has an approximate value of 0.0307, but the true value of w is 
See ca 
9+ 11+°°" 


3. If 0.54931 = x < 9.011, then use the identity tanh (x) = 1 — <7 The 


computation for this case uses the real exponential subprogram (IHCSEXP). 
4, If x = 9.011, then tanh (x) = 1. 
5. If x S —0.54931, then use the identity tanh (x) = —tanh (—~). 


Effect of an Argument Error 


E ~ (1 — tanh? x) A, and e . For small values of x, « ~ 8, and as the 


2A 
sinh (2x) 
value of x increases, the effect of § upon e diminishes. 


IHCSTNCT Subprogram (TAN and COTAN) 
Algorithm 


1. Divide |x| by-jand separate the result into the integer part (q) and the 


fraction part (17). Then, 





x| = 4 + r), 


2. Obtain the reduced argument (w) as follows: 
if q is even, then w = 1, 
if gq is odd, then w = 1 — 1. 
The range of the reduced argument is 0 = w = 1. 


3. Let go = q mod 4. 
Then, for qo = 0, tan |x| = tan ( 4 . w) and cot |x| = cot ( ri . w) 
. w) 














4 


qo = 1, tan |x| = cot ( . w) and cot |x| = tan ( z 








Jo = 2, tan |x| = —cot ( ji ° w) and cot |x| = —tan ( i ° w) 


go = 3, tan |x| = —tan ( r ° w) and cot |x| = —cot (+ . w) 











4. The values of tan ( ri ° w) and cot ( r . w) are computed as the ratio of 


two polynomials. 
7. VL wt P(w?) TS wp) x Oi) 
fore a) oa) 
where P(w?) = 212.58037 — 12.559912w? 
Q(w?) = 270.665736 — 71.645273w? + wt 
This approximation is obtained by economizing the continued fraction 
tan(z)_ 4) _ #2 2 zt 
z 3- 5-— 7-""° 
in the following way: 


2 2 2 
Write: tan(z) Zz Zz z 








tan 





eS ee 


(3+d,)— (5+ d2.)— (7+ ds) 
and determine values for d,, dz, and d3 so that the right-hand expression 
gives the exact answers for z2 = 0.19, 0.432, and 0.594. Then the maximum 





relative error of this formula over the range 0 Sz 7 1.74 + 10-8. 








Change the variable from z to w = ao, z and rewrite the formula to ob- 
T 
tain P(w?) and QO(w?). 
5. If x < 0, then tan |x|, and cot(x) = —cot |x]. 
Effect of an Argument Error 
2 ; She 
aoa and « ~ sin( x) for tan(x). Therefore, near the singularities 


1 : ets 
x= (: +z) where k is an integer, no error control can be maintained. 


This is also true for cotan(x) for x near kz, where k is an integer. 
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Appendix B. Performance Statistics 


Appendix B contains accuracy and timing statistics for 
the explicitly called mathematical subprograms. These 
statistics are presented in Table 12 and are arranged 
in alphabetical order, according to the entry names. 
The following column headings are used in Table 12: 


Entry Name: This column gives the entry name that 
must be used to call the subprogram. 


Argument Range: This column gives the argument 
range used to obtain the accuracy figures. For each 
function, accuracy figures are given for one or more 
representative segments within the valid argument 
range. In each case, the figures given are the most 
meaningful to the function and range under consid- 
eration. 


The maximum relative error and standard deviation 
of the relative error are generally useful and revealing 
statistics; however, they are useless for the range of a 
function where its value becomes 0, because the slight- 
est error in the argument can cause an unpredictable 
fluctuation in the magnitude of the answer. When a 
small argument error would have this effect, the maxi- 
mum absolute error and standard deviation of the 
absolute error are given for the range. For example, 
absolute error is given for sin(x) for values of x near 7. 


Sample: This column indicates the type of sample 
used for the accuracy figures. The type of sample de- 
pends upon the function and range under consid- 
eration. The statistics may be based either upon an 
exponentially (E) distributed argument sample or a 
uniformly (U) distributed argument sample. 


Accuracy Figures: This column gives accuracy fig- 
ures for one or more representative segments within 
the valid argument range. The accuracy figures sup- 
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plied are based upon the assumption that the argu- 
ments are’ perfect (i.e., without error and, therefore, 
having no error propagation effect upon the answers). 
The only error in the answers are those introduced by 
the subprograms. Appendix A contains a description 
of some of the symbols used in this appendix; the 
following additional symbols are used in the presenta- 
tion of accuracy figures: 


f(x) — g(x) 
f(x) 


The maximum 
relative error 
produced 
during testing. 


The standard 

3 deviation (root- 
mean-square ) of 
the relative 
error. 


Max 


M (ce) 








] f(x) — g(%) 
a (e) wet £(x;) 





The maximum 
absolute error 
produced 

during testing. 


The standard 
deviation (root- 
mean-square ) of 
the absolute 
error. 


M (E) 


Max | f(x) — g(x) | 


F(x) — g(a) 


o(B) = Ia; 








In the formulas for the standard deviation, N repre- 
sents the total number of arguments in the sample; 
iis a subscript that varies from 1 to N. 

Average Speed: This column gives the timing statis- 
tics. These statistics represent the average speed in 
microseconds for the various System/360 models. Sta- 
tistics are supplied for Models 30, 40, 50, 65, and 75. 
Statistics for the Model 75 are based upon two-way 
interleaving. 


sonsieig souvullojieg ‘g xipueddy 


Table 12. Performance Statistics 





ce daa | ee 
Name Range E/U M (e) ete M (EB) o (E) 30 40} 50} 65 75 
(See Note8) 
ALGAMA | 0< x05 U 1.09 x 10-6 3.35 X 10-7 10500} 2800] 865] 221 131 
| foscxca fw | dYSS*dS eto? | ata 10-7 | 10500] 2520) 664] 295 155 | 
3=x<8 U 1.21 x 10-6 2.86 x 10-7 12100] 3250|1020| 259 151 


= 5] 
> Q 
Z fe) 

172) 


ala 


saxcie | _U | 1asxi0-* | saoxio-7 | ——S«dT—SSSCSCS~S~C~CS 
jesx<sw [Uv [roixao-e | a0axio-? [|_| | revo] aorol erz|isa of 
055x515 | u {| | ss 346 x 10-7 =| 8.62 x 10-8 91.9 52.3 
r<05x>15 | E | aaxw-* | ia KIO | 261[ O19 52. 

PP 64 x 10-7 4.78 X 10-8 388| 98.1 56.1 
ALOGIO 66.1561 
P-igsaet 


Gg 


ALOG 






te 


1.80 x 10-7 | 3.29 x 10-6 5340| 1460] 451|118 70.6 
SIN U 8.56 x 10-7 | 2.38 x 10-7 5270| 1450] 445/114 68.4 


The full range Note 7 9.75 x 10-7 4.54 x 10-7 For ATAN in FORTRAN IV (E) 


3602 913 255 68.8 40.6 
ATAN2 | The full range 9.75 x 10-7 | 4.54 x 10-7 
















5.9] 
9.7 


Additional time for ATAN in 
CABS The full range 1.87 x 10-6 | 7.65 x 10-7 


FORTRAN IV 
CCOos |x: | = 10, |xo| = 20 U 1.79 x 10-6 7.21 X 10-7 
See Note 2 
<= U 


104 24 
4874| 1288 1031 
CDABS | The full range 3.32 x 10-15 | 5.16 x 10-16 
CDCOS | |x| < 10, |xe| <1 5.16 X 10-15 | 3.42 x 10-16 
See Note 3 
CDEXP | |x| = 1, |x| S - 4.04 X 10-16 | 1.39 x 10-16 


5194] 1421] 396/108 68.4 
15783 | 4433| 132 203 
|x| = 20, |xe| S 20 U 3.63 x 10-15 | 1.29 x 10—16 
CDLOG | The full range 8.73 X 10-15 | 6.38 x 10-17 


48343 |11705 | 2335 
CDSIN |xa| = 10, |x.| =1 U 3.72 X 10-15 3.49 X 10—16 
See Note 4 


50 89.0 
301 


41737 |10144|2048| 456 260 


7 
48250 |11668 |2322|503 298 
1 


27951 | 5906] 1282 18 


13731 | 3888;1166; 291 17 


13899 | 3930/1180] 294 17 
15504 | 4246)1261| 338 21 


— 
i 
oO 
bo 
_ 
(ey) 
So 
lop) 
_— 
> 
ive) 
(ee) co) 


H ) 
Oo 
Pec 





cl 
bw. 


CDSQRT | The full range Note 1 9.86 x 10-16 | 1.91 x 10-16 


CEXP | |u| $170, fu 37 boo 11s xi0-6 | 9.34% 10-7 


U 
[xa] = 170, 
F< 
CLOG The full range 
except (1 + 0i 


7 
8 
6 


OF 







Table 12. Performance Statistics (Continued ) 


Accuracy Figures Average Speed ( Microseconds ) 
E/U 


(See Note 8) 
7 — 












a 1.47x 10-7 | 5.48 x 10-8 3934 | 1047 742| 44.0 
Ld 1.42 X 10-7 5.67 xX 10-8 3990 | 1061 | 303] 75. 


A 
~ 
a 


aN 
» 
oO 


~I ~I 
on on 
i m| bo 








o<haw | 0 


1.31 x 10-6 


oO 
co 
bo 


3.40 x 10—6 6110 | 1810 | 570] 145 


4420 | 1180 | 341 | 86. 


eae ae 1.35 X 10-7 5.61 X 10—8 3990 | 1061 | 303] 75. 


~] 
Shi 
fon) 
Oo 


Le 15690 | 4397 [1316 | 331 | 200 


aN 
fon) 
— 
fom) 
jou 
po 
bo 
=) 
KR 
— 
oe) oO 
oO ce) 
_ Go 
Ut wt 
or Ol 
Oo Oo 


co 
io) 
jt 
91) 
OL 
i) 


4 
ae = 1.13 x 10-5 6.03 X 10-7 
2 < | 0 See Note 5 
10 < |x| S 100 1.67 X 10-5 6.67 X 10-7 
See Note 5 
CSIN |x:] S 10, [xe] S 1 1.97 X 10-& | 7.09 x 10-7 
See Note 6 


CSQRT The full range Note 1 1.61 x 10-6 4.58 X 10-7 
2.72 X 10-16 | 9.35 x 10-17 


For DATAN in FORTRAN IV(E) 
19056 4000 715 153 83.6 
Additional time for DATAN in 
FORTRAN IV 


U 
ax 
COTAN |x] = = 129x 10-8 | 3.68 x 10-7 a 
wT <7 3.80 x 10-4 7.70 X 10-4 
sss | 8 [ewes | | 
: hod 
: aaa 
: booed 
U 










The full range Note 7 2.08 X 10-16 | 6.64 x 10-17 a 


The full range 2.08 x 10-16 | 6.64 x 10-11 
1.79 x 10-16 | 6.40 x 10-17 


[peered 





DATAN2 


-103x<0 
r<x=l0 
hlo<jx}=100 =|) 6ChUU Ts 865 10-15 | 1.01 x 10-15 13133 | 3146 | 605| 132. | 


U 

U 

U 
i : = 18300 | 4260 | 898 

DCOTAN | |x| <7 . 3.46 X 10-16 | 8.38 x 10-17 

U 1.72 X 10-18 | 5.00 x 10-15 
See Note 5 
= U 5.33 X 10-13 | 1.09 x 10-14 






10 < |x| = 100 8.61 X 10-13 | 4.61 x 10-14 
See Note 5 
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Table 12. Performance Statistics (Continued ) 


a | heme eat Average Speed ( Microseconds) 


Name Range E/U M (ce) pts) M (E) o (E) 30 40} 50] 651 75 
(See Note 8) 
DERF [x] S 1.317 1.70 x 10-16 | 2.71 x 10-17 18400! 4610! 8791 190 | 109 


1.317 < |x| S 2.04 
2.04 < |x| << 6.092 
-6<x<0 
0=x=1.317 
1.317 << x S 2.04 
204<x<4 
4=x<133 
DEXP xj S1 

1 < |x| S 20 

20 < |x| = 170 
O0<x<l 
b=a7=2 


2<x84 


~[aerxa0= Paarxao |i 200 60 250 aa 
1.70 X 10-17 | 8.03 x 10-18 | _ 45200|10900|2080| 449 | 259 
[iss «10-16 | oa eao- [+d SSS 8800 
ase x10" | 762 x10-"# [| ——«P——~ 0 260 
02 x 10- ces 
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aati | 7a9xao- [dT SS—*S af 0 

asrx 108 | 360 a0 [iY ia 2 

esx [essxio— | «| _————~d—a as | aor] 
et! 
een 
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OoO;1o 
aja 
— | 
Ww) vw 
GO | GO 
~1j 1 
oul Ot 
or] ot 


aq 
—_ 
w 
Go 
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DGAMMA 


218x 10-16 | 7.93x10-77 | 30700 | 7500 {140 

3.12x 10-17 | 845x10-18 | 28200| 7050 

169x 10-15 | 460x10-79 | | sf 30100] 7520]14 
3 


o 


~] 
oe 
o 
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«J 
Ot 


bt | ee 
aS 
oO;o 
wyhe 
bet | CO 
BE 
a 
~wloa 
EG 
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So 
oO 
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oO 
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jad 
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=] 


4<x<8 2.85 X 10—15 9.46 X 10-16 3900 | 8470 |1610} 352 | 205 
ean 
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05 <x<8 Pf B10 ae TOS 00 fa 
pax<8 [saxo [iisxaoe | 
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DLOG10 
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7 
x <0.5,x>15 6.14 x 10-16 | 9.96 x 10-17 i ee 17147 | 4044 170 | 92.0 
DSIN kis 4.08 x 10-16 | 4.85 x 10-17 | 9.10x 10-17 | 9.17 x 10-27 13145 | 3148 133 | 75.5 


Tr 


| 

A 

a 

lA 

S 
cf 


ae 1.64 x 10-16 6.35 X 10-17 13145 | 3148 


ecg wee ol 2.69 X 10-15 1.03 X 10-15 13145 | 3148 


10 < |x|S 100 


DSINH | |x| < 0.34657 
0.34657 < |x| <5 
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Sample Accuracy Figures Average Speed { Microseconds ) 


Table 12. Performance Statistics (Continued ) 
E/U M (e) glee M (E) o (E) 30 40 | 50] 65 1 75 
(See Note 8) 


Entry 
Name ) 
The full range E 108x 10-16] 2i7x10-7 | = ss—i‘iLss—s~sdY~Sseit | 20084 | 355] 85.3] 49.2 


Argument 
Range 







DSQRT 
aap <7 U 1.67 X 10-12 | 3.69 x 10-14 15700 | 3660 
4 = 9 See Note 5 
1.57 x 10-13 | 4.51 x 10-15 15600 | 3640 
See Note 5 
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S 
oo 


Ee 


| 8 
oO 
oo 


— 
Ot 
oO 
[ee] 
fon) 
—_ 


v 

me < 

9 < |x| S10 

10< |x| = 100 U 3.79 X 10-12 | 9.50 x 10-14 
See Note 5 


0.54931 < |x| $5 


[x] = 170 


ERF |x| S 1.317 
1.317 < |x| S 2.04 
2.04 < |x| S 3.9192 


ERFC -38<x<0 
0Sx51317 
1317 <x =2.04 
: : 


DTANH | |x| < 0.54931 
199 x 10-16 | 2.54 x 10-17 


4.65 x 10-7 1.28 x 10—-T é : 
4.69 x 10-7 1.17 X 10-7 4183 | 1250 | 387 | 94. 


GS 
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9.26 x 10-7 1.43 x 10-7 


2 
pss 
oO 
— 
co) 
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OQ 
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gt 
(vy) 
co 
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ig 


374 | 92.0 ae 


0| 58.0 
210 | 126 


eorxio [eats | 
eax [aorxi® | 
seo | sax | 
pueexio* [aisxioF | 
pisoxi0 [300x107 | 
= ans 

36x 10- ae: 
Siac! 

ae 

aaa 


Pa 

~~ 

—_ 

oO 
E 

bo 

_ 

a] 
BA 

=] 


ilsee 

[aoexi0-* [eax [ «it ~*~ 0 | 
aor 
; anaes) 


8 
9950 | 2780 | 835 
1.52 x 10-5 8.45 x 10—6 1 
9.86 x 10-7 3.45 X 10-7 
1.00 x 10-7 3.74 X 10-8 


=| 8.4 x 10- 77a0 
2299x10-5 | 767x10-@ | ss 12000 | 3320 |1020 | 263 | 155 
436x10-5 | 145x10-5 [ | ss 12000 | 3320 f1020] 263 | 155 
159x 10-8 | 202x10-7 | 131x 10-7 | 5.55 x 10-8 3876 | 1036 | 298| 74.2 


1.41 xX 10-7 5.53 X 10-8 3989 | 1064 | 307} 76.4} 45. 


2.00 x 10-16 | 4.45 x 10-17 


as < 
7) < |x| = 10 


10 < |x| S 100 
SQRT The full range 
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Table 12. Performance Statistics (Continued ) 
Accuracy Figures 


Average Speed (Microseconds ) 
30 40 50 i 75 


6 
(See Note 8) 


3.22 xX 10-7 





1.56 x 10-6 


7 U 6.58 X 10-5 
2 See Note 5 


wW<|x}<100 | UU 3.35 X 10-5 
See Note 5 
TANH |x| S 0.54931 8.12 x 10-7 | 1.66 x 10-7 
0.54931 < |x| $5 5.74 x 10-7 | 7.53 x 10-8 


1.67 x 10-6 4500 | 1184 | 338) 85.3) 52.8 


Note 4: The maximum relative error cited for the 
cpsin function is based upon a set of 1500 random 
arguments within the range. In the immediate prox- 
imity of the points nmr + Oi (where n = +1, +2, 
...,) the relative error can be quite high although 
the absolute error is small. 


Notes to Table 12 


These notes are associated with Table 12 and contain 
more detailed information about samples and relative 
errors for certain functions. 


Note 1: The distribution of sample arguments upon 
which these statistics are based is exponential radially 
and is uniform around the origin. 


Note 5: The figures cited as the maximum relative 
errors are those encountered in a sample of 2500 ran- 
dom arguments within the respective ranges. See the 
appropriate section in Appendix A for a description 
of the behavior of errors when the argument is near a 
singularity or a zero of the function. 


Note 2: The maximum relative error cited for the 
ccos function is based upon a set of 2000 random argu- 
ments within the range. In the immediate proximity of 


Note 6: The maximum relative error cited for the 
csIN function is based upon a set of 2000 randon argu- 
ments within the range. In the immediate proximity 
of the points nz + 0i (wheren = +1, +2,...,) the 
relative error can be quite high although the absolute 
error is small. 


the points (» + +) a + Oi (where n= 0, £1, +2, 


...,) the relative error can be quite high, although 
the absolute error is small. 


Note 3: The maximum relative error cited for the 
cpcos function is based upon a set of 1500 random 


arguments within the range. In the immediate prox- Note 7: The sample arguments were tangents of 


numbers uniformly distributed between — a and + oe 


SOTJSI}FVIG VoULUTIOJAeg ‘gq xIpueddy 


cP 


imity of the points (» + +) ax + Oi (where n = 0, 


#1, +2,...,) the relative error can be quite high 
although the absolute error is small. 


Note 8: The statistics for the Model 75 are based 
upon two-way interleaving. 


Appendix C. Interruption and Error Procedures 


Appendix C contains descriptions of the procedures 
followed when the execution of a load module is dis- 
continued. Execution may be discontinued due to one 
of two reasons: an interruption or an error. After an 
interruption is processed, execution of this load module 
continues; after an error is processed, execution of this 
load module is terminated. The following text explains 
the procedure used to handle each case. 


Interruption Procedures 


An interruption is a computer-originated break in the 
flow of processing. An interruption occurs when an 
invalid arithmetic operation is attempted or when a 
result which cannot be fully contained in a floating- 
point register is returned. These interruptions are due 
specifically to one of the three following exceptions: 

1. Exponent-overflow exception. 

2. Exponent-underflow exception. 

3. Divide exception. 

An exponent-overflow exception is recognized when 
the result of a floating-point addition, subtraction, 
multiplication, or division.is either greater than or 
equal to 16®* (approximately 7.2 X 107). An expo- 
nent-underflow exception is recognized when the result 
of a floating-point addition, subtraction, multiplication, 
or division is less than 16—-® (approximately 5.4 X 
10-79). A divide exception is recognized when division 
by zero is attempted. 

For all three exceptions, an indicator is set and a 
message is written in the data set associated with 
system output. After the interruption is handled, the 
execution of this load module continues from the point 
of interruption. 

The status of an indicator corresponds to the excep- 
tion that caused the interruption. The indicators can 
be tested with the subprograms described in the sec- 
tion, “Service Subprograms.” 

The message contains the old program status word 
(esw) and has the format shown in Figure 1. One of 


IHC2101 


Figure 1. Program Interrupt Message Format 
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PROGRAM INTERRUPT--OLD PSW IS xxxxxxx 


four characters (9, C, D, or F) appears in the old psw 
and indicates the cause of the interruption, as follows: 
9 — fixed-point divide exception 
C — exponent-overflow exception 
D — exponent-underflow exception 
F — floating-point divide exception 
More information on the psw can be found in the pub- 
lication IBM System/360: Principles of Operation. 


Error Procedures 


During execution, the mathematical subprograms as- 
sume that the argument(s) is the correct type. No 
checking is done for erroneous arguments (i.e., the 
wrong type, invalid characters, the wrong length, etc. ); 
therefore, a computation performed with an erroneous 
argument has an unpredictable result. However, the 
nature of some mathematical functions requires that 
the input be within a certain range. For example, the 
square root of a negative number is not permitted. 
If the argument is not within the valid range given in 
Tables 2 through 6, an error message is written in the 
data set associated with system output. The execution 
of this load module is terminated and control is re- 
turned to the operating system. 

The error message that is issued has the format 
1Hcxxx1. The xxx is a numeric code that identifies the 
error detected. The following text lists the error mes- 
sages in numeric order, explains the error, and indi- 
cates what action the system takes. In the following 
explanations, x represents the argument supplied by 
the programmer. 


1HC216I 


Explanation: In the micrstir subprogram, a value of 
i that is not 0, 1, 2, 3, or 4 is an error. 

System Action: Execution of this load module is 
terminated. 


XXXXXKKXX 





1HC2411 
Explanation: In the racrixp1 subprogram, a base 
number of zero and an exponent = 0 is an error. 
System Action: Execution of this load module is 
terminated. 


1HC242! 
Explanation: In the rHcrrxP1 subprogram, a base 
number of zero and an exponent = 0 is an error. 
System Action: Execution of this load module is 
terminated. 


1HC2431 
Explanation: In the micrpxp1 subprogram, a base 
number of zero and an exponent <= 0 is an error. 
System Action: Execution of this load module is 
terminated. 


1HC244I 
Explanation: In the tHcFRxPR subprogram, a base 
number of zero and an exponent = 0 is an error. 
System Action: Execution of this load module is 
terminated. 


1HC245! 
Explanation: In the micrpxpp subprogram, a base 
number of zero and an exponent = 0 is an error. 
System Action: Execution of this load module is 
terminated. 


1HC246! 
Explanation: In the tHcrcxr1 subprogram, a base 
number of zero and an exponent = 0 is an error. 
System Action: Execution of this load module is 
terminated. 


1HC247I1 
Explanation: In the racrcpx1 subprogram, a base 
number of zero and an exponent = 0 is an error. 
System Action: Execution of this load module is 
terminated. 


JHC2511 
Explanation: In the macssert subprogram, a value 
of x < Ois an error. 
System Action: Execution of this load module is 
terminated. 


1HC2521 


Explanation: In the micsExp subprogram, a value of 
x > 174.673 is an error. 

System Action: Execution of this load module is 
terminated. 


IHC2531 


Explanation: In the t9csLoc subprogram, a value of 
x = 0 is an error. Because this subprogram is also 


called by an exponentiation subprogram, this mes- 
sage also indicates that an attempt has been made 
to raise a negative real base to a power. 

System Action: Execution of this load module is 
terminated, 


1HC254! 
Explanation: In the 1Hcsscn subprogram, a value of 
|x|= 28 + @ is an error. 
System Action: Execution of this load module is | 
terminated. 


1HC255! 
Explanation: In the tHcsatne subprogram when 
entry name ATAN? is used, a value of x1 = x2 = 0 is 
an error. 
System Action: Execution of this load module is 
terminated. 


1HC256! 
Explanation: In the t1csscnH subprogram, a value 
of |x| = 174.673 is an error. 
System Action: Execution of this load module is 
terminated. 


1HC2571 
Explanation: In the macsascn subprogram, a value 
of |x| > 1 is an error. 
System Action: Execution of this load module is 
terminated. 


1HC258!1 
Explanation: In the ticstNcr subprogram, a value of 
|x| = 218 « + is an error. 
System Action: Execution of this load module is 
terminated. 


1HC2591 
Explanation: In the racstNcr subprogram, a value of 
a 3 
x too close to one of the singularities ( + a ae va 


for the tangent; +7, +27,... for the cotangent) is 
an error. 
System Action: Execution of this load module is 


terminated. 


1HC2611 , 
Explanation: In the tHcLsgrt subprogram, a value 
of x < 0 is an error. 

System Action: Execution of this load module is 
terminated. 


1HC2621 


Explanation: In the tHcLExp subprogram, a value of 
x > 174.673 is an error. 

System Action: Execution of this load module is 
terminated. 
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1HC2631 


Explanation: In the 14cLLoc subprogram, a value of 
x <= 0 is an error. Because this subprogram is also 
called by an exponentiation subprogram, this mes- 
sage also indicates that an attempt has been made 
to raise a negative real number to a power. 

System Action: Execution of this load module is 
terminated. 


1HC2641 


Explanation: In the 19cLscn subprogram, a value of 
|x| 2 25° + 7 is an error. 

System Action: Execution of this load module is 
terminated. 


1HC2651 
Explanation: In the 1HcLATN2 subprogram whe. 
entry name DATAN? is used, a value of x; = x2 = 0 
is an error. 
System Action: Execution of this load module is 
terminated. 


IHC266I 


Explanation: In the 1cLscNH subprogram, a value 
of |x| => 174.673 is an error. 

System Action: Execution of this load module is 
terminated. 


1HC2671 


Explanation: In the tcLAscn subprogram, a value 
of |x| > 1 is an error. 

System Action: Execution of this load module is 
terminated. 


1HC2681 
Explanation: In the ncLtNncr subprogram, a value of 
|x| = 25° + 7 is an error. 
System Action: Execution of this load module is 
terminated. 


1HC2691 
Explanation: In the 1ucLtNct subprogram, a value of 
x too close to one of the singularities ( + ae +, . 


for the tangent; 7, +27, ... for the cotangent) is 
an error. 
System Action: Execution of this load module is 


terminated. 


IHC2711 


Explanation: In the miccsExp subprogram, a value 
of x, > 174.673 is an error. 

System Action: Execution of this load module is 
terminated. 
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1HC2721 
Explanation: In the miccsExp subprogram, a value 
of |x| 2 218 » 7 is an error. 
System Action: Execution of this load module is 
terminated. 


1HC273! 
Explanation: In the miccsLoc subprogram, a value 
of x; = x2 = Ois an error. 
System Action: Execution of this load module is 
terminated. 


1HC274I 
Explanation: In the raccsscn subprogram, a value 
of |2x1| => 218 » 7 is an error. 
System Action: Execution of this load module is 
terminated. 


1HC275I 
Explanation: In the taccsscn subprogram, a value 
of |x2|> 174.673 is an error. 
System Action: Execution of this load module is 
terminated. 


IHC2811! 
Explanation: In the miccLExpe subprogram, a value 
of x, > 174.673 is an error. 
System Action: Execution of this load module is 
terminated. 


1HC282I 


Explanation: In the taccLexp subprogram, a value 
of |x| 2 950 » is an error. 

System Action: Execution of this load module is 
terminated. 


1HC283! 
Explanation: In the ticcLLoc subprogram, a value 
of x; = x, = Ois an error. 
System Action: Execution of this load module is 
terminated. 


1HC284! 
Explanation: In the miccLscn subprogram, a value 
of |21| = 250 « + is an error. 
System Action: Execution of this load module is 
terminated. 


1HC285! 
Explanation: In the tHccuscn subprogram, a value 
of |x2| > 174.673 is an error. 
System Action: Execution of this load module is 
terminated. 


IHC290I 


Explanation: In the 19csGAMA subprogram for the 
gamma function, a value of x = 2-5? or x = 57.5744 
is an error. 


System Action: Execution of this load module is | 


terminated. 


tHC2911 


Explanation: In the rucsGAMA subprogram for the 
log-gamma function, a value of x = 0 or x = 4.2937 
* 1078 is an error. 

System Action: Execution of this load module is 
terminated. 


IHC3001 


Explanation: In the mtcLGAMA subprogram for the 
gamma function, a value of x SS 2-252 or x = 57.5744 
is an error. 
System Action: Execution of this load module is 
terminated. 


IHC3011 


Explanation: In the t9cLGAMA subprogram for the 
log-gamma function, a value of x = 0 or x = 4.2937 
¢ 1078 is an error. 

System Action: Execution of this load module is 
terminated. 
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Appendix D. Storage Estimates 


Appendix D contains decimal storage estimates (in 
bytes) for the library subprograms. The estimate given 
does not include any additional library subprograms 
Or FORTRAN execution-time routines that the subpro- 
gram needs during execution. The names of any addi- 
tional library subprograms needed are given in Tables 
13 and 14 in the column headed “Additional Subpro- 
grams.” 


Some library subprograms also require execution- 
time routines for input/output, interruption, and error 
procedures. The IHCFIosH routine performs input/ 
output procedures for both FoRTRAN Iv (E) and FORTRAN 
1v. (This routine refers to a table (1HcUaTBL) for infor- 
mation about the input/output devices used during 
execution. ) The 1mcFCcoME routine performs interrup- 
tion and error procedures for FORTRAN Iv (E) library 
subprograms; the IHCFCOMH-IHCFCVTH routines perform 
these procedures for FoRTRAN Iv library subprograms. 
If an installation has both compilers, the micrcoMH- 
IHCFCVTH routines are used. Tables 13 and 14 indicate 
which library subprograms require these execution- 
time routines. 


In addition, several other execution-time routines 
may be needed to resolve external references in a 
FORTRAN Iv object module. 


1. If a source module specifies direct-access input/ 
output operations, the compiler generates a call to 
the THCDIOSE routine. 


2. At the point that errors are encountered during 
compilation, the compiler generates a call to an 
error routine (IHCIBERR for FORTRAN Iv (E) and 
IHCIBERH for FORTRAN Iv). If execution of the load 
module is attempted, the error routine is called, 
a message is issued, and execution is terminated. 


3. If a FORTRAN Iv (£) source module contains a com- 
puted co To, the compiler generates a call to the 
IHCCGOTO routine. 


4, If a FORTRAN IV source module contains any input/ 
output operations that refer to a NAMELIST name, 
compiler generates a call to the 1#CNAMEL routine. 


5. If a FORTRAN IV source module uses the debug 
facility, the compiler generates a call to the ncpBUG 
routine. 
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Table 13. Mathematical Subprogram Storage Estimates 


Subprogram 
Name 


IHCCLABS 
THCCLAS 
IHCCLEXP 
IHCCLLOG 


IHCCLSCN 
IHCCLSQT 
IHCCSABS 
IHCCSAS 
IHCCSEXP 
IHCCSLOG 


IHCCSSCN 
IHCCSSQT 
IHCFAINT 
IHCFCDXI 
THCFCXPI 
JHCFDXPD 
IHCFDXPI 
IHCFIFIX 
JHCFIXPI 
IHCMFAXD 
IHCFMAXI 
IHCFMAXR 
ITHCFMODR 
IHCFMODI 
IHCFRXPI 
IHCFRXPR 
IHCLASCN 
IHCLATAN 
IHCLATN2 
IHCLERF 
IHCLEXP 
IHCLGAMA 
IHCLLOG 
IHCLSCN 
IHCLSCNH 
IHCLSQRT 
IHCLTANH 
IHCLTNCT 
IHCSASCN 
IHCSATAN 
IHCSATN2 
IHCSERF 
IHCSEXP 
IHCSGAMA 
IHCSLOG 
IHCSSCN 
IHCSSCNH 
THCSSQRT 
IHCSTANH 
IHCSTNCT 


Uses 
IHCFIOSH 
and 
IHCFCOME 


Additional 
Subprograms 


IHCLSQRT 


IHCLEXP, IHCLSCN 
IHCCLABS, IHCLSORT, 
IHCLLOG, IHCLATN2 
IHCLEXP, IHCLSCN 
IHCLSORT 

IHCSSORT 


IHCSEXP, IHCSSCN 
IHCSABS, IHCSSQRT, 
IHCSLOG, IHCSATN2 
IHCSEXP, IHCSSCN 


THCSSQRT 


THCCLAS 
IHCCSAS 
IHCLLOG, IHCLEXP 


IHCSLOG, IHCSEXP 
IHCLSQRT 


IHCLEXP 


IHCLLOG, IHCLEXP 


IHCLEXP 
IHCLEXP 
IHCSSQRT 


IHCSEXP 


IHCSLOG, IHCSEXP 


IHCSEXP 


IHCSEXP 





If the programmer has not made allowances for the 
storage required by any of these additional routines 
(see Table 15), the amount of available storage may 
be exceeded and execution cannot begin. The pro- 
grammer must add the estimates for all subprograms 
and routines needed to determine the amount of stor- 
age required. 


Table 14, Service Subprogram Storage Estimates 


Uses IHCFIOSH and 
IHCFCOME or 
IHCFCOMH-IHCFCVTH 












Decimal 
Estimate 


Subprogram 
Name 


IHCFDVCH 
IHCFDUMP 
IHCFEXIT 
THCFOVER 
IHCFSLIT 








Table 15. Execution-Time Routine Storage Estimates 


Routine Decimal | 
Name Estimate Used By 
























































IHCCGOTO 60 FORTRAN IV (E) 
IHCDBUG 2,600 FORTRAN IV 
IHCDIOSE 2,500 Both 
IHCFCOME 5,500 FORTRAN IV (EB) 
IHCFCOMH 4,050 FORTRAN IV 
IHCFCVTH 4,090 FORTRAN IV 
IHCFIOSH 3,800 + IHCUATBL Both 
(See Note) 

IHCIBERH 210 FORTRAN IV 
IHCIBERR 260 FORTRAN IV (E) 
IHCNAMEL 2,250 FORTRAN IV 





Note: The number of bytes in table IHCUATBL may be 
computed by the formula 

12n + 8 
where n is the number of data set reference numbers 
requested during system generation. 
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Appendix E. 


The mathematical and service subprograms in the 
FORTRAN Iv library are available to the assembler lan- 
guage programmer. The following text explains the 
method of calling a library subprogram in an assembler 
language program, and then gives additional informa- 
tion necessary to use each type of subprogram. (The 
assembler language programmer should also be fami- 
liar with the information contained in Appendix D.) 


Calling Sequences 


To call either type of library subprogram, the assem- 
bler language programmer supplies an entry name, 
an argument list, and an area used by the subprogram 
to store information (i.e., a save area). The following 
conventions must be observed when calling a library 
subprogram in an assembler language program: 

1. The address of the entry name must be in general 
register 15. 

2. The address of the point of return to the calling 
program must be in general register 14. 

3. The address of the argument list must be in general 

register 1. 

4, The argument list must be assembled on a full-word 
boundary; it consists of one 4-byte address constant 
for each argument. The last argument must have 
a 1 in its high order bit. 

5. The address of the save area must be in general 
register 13. 

6. The save area must be assembled on a full-word 
boundary. Although the minimum size of the save 
area depends upon the subprogram, the program- 
mer is advised to use a save area of 18 full-words 
for all library subprograms. The minimum save 
area sizes are given in Tables 2 through 6 for the 
mathematical subprograms, and in Table 16 for the 
service subprograms. 

7. If the information in a floating-point register is to 
be retained, the programmer must save and restore 
the contents of the register. The subprograms that 
make use of the floating-point registers contain no 
provisions for saving the information. 

8. If a main program in assembler language contains 
any calls to those library subprograms that use the 
FORTRAN execution-time routines (see Appendix D), 
the following instructions must be included before 
the call to the subprogram is issued: 

L 15, = V(IBCOM#) 

BAL 14,64(15) 

These instructions cause the initialization of return 
coding and the interruption exceptions described in 
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Assembler Language Information 


Appendix C. If these instructions are omitted, the 
occurrence of an interruption or an error causes un- 
predictable termination of the execution of this load 
_ module. 
Note: In an assembler language program, a decimal- 
divide exception may occur. This causes the char- 
acter B to appear in the program interruption 
message described in Appendix C. 


The assembler language programmer may use sev- 
eral methods to call a ForTRAN library subprogram: the 
appropriate macro-instructions described in the pub- 
lication, IBM Operating System/360: Control Program 
Services or the general assembler language calling 
sequence (given in Figure 2). If the macro-instruc- 
tions are used, the address of the save area must be 
placed in general register 13 before using a macro- 
instruction to give control to the subprogram. For 
example, if the square root of the value in AMNT is to 
be computed and sAvE is the address of the save area, 
the following statements could be included in an 
assembler language program to call the maicssgrt sub- 
program: 


L 15,=V(IBCOM#) 
BAL 14,64(15) 

LA " 13.SAVE 

CALL SORT,( AMNT),VL 


SAVE DS ' 18F 

If the general assembler language calling sequence 
shown in Figure 2 is used, the programmer must en- 
sure that all of the conventions discussed previously 
are followed. For example, to call the macssgrr sub- 
program to compute the square root of the number in 
AMNT, the following statements would be included in 
the source program: 


LA 13,SAVE 
LA LARG 
L 15, ENTRY 
BALR 14.15 
ENTRY DC V(SORT) 
SAVE DS  18F 
ARG DC X‘80" 
DC AL3(AMNT) 


LA 13, area 
LA 1, arglist 
L 15, entry 
BALR 14, 15 
NOP X’id’ 


* * * * 


DC -V (entry name) 
or 


DC A (entry name) 
* * * 


DS xxF 


* * * * 


CNOP 
arglist ‘DC X’80’ 
DC AL3 (arg) 
or for more than one argument: 
arglist DC A (arg; ) 
DC A (argz) 


DC X’80’ 
DC AL3 (argn) 


Figure 2, General Assembler Language Calling Sequence 








When the load module is executed the rucssert sub- 
program is called to compute the square root of the 
number in AMNT; the result is stored in floating-point 
register 0. The binary calling sequence identifier is 
not used. 


Mathematical Subprograms 

The assembler language programmer supplies one or 
more arguments for each mathematical subprogram. 
The arguments may be either integer values or nor- 
malized floating-point real or complex values. 

An integer argument is four bytes in length and 
starts on a full-word boundary. A real argument is 
either four or eight bytes in length. The four-byte 
argument starts on a full-word boundary. The eight- 
byte argument starts on a double-word boundary and 
occupies two adjacent words. The first word contains 
the most significant digits. This word is also the ad- 
dress of the entire argument; the second word con- 
tains the least significant digits. 

A complex argument is either eight or sixteen bytes 
in length and starts on a double-word boundary. The 
first half of the argument contains the real part of the 
complex argument; the second half contains the 
imaginary part. The address of the real part of the 
argument is the address of the entire argument. 

Each mathematical subprogram returns a single 
answer. This answer is either an integer value or a 


General register 13 contains the address of the save area. 

General register 1 contains the address of the argument list. 

General register 15 contains the address of the subprogram. 

General register 14 contains the address of the point of return to the 
calling program. 

This statement is optional. The id represents the binary calling sequence 
identifier. This number is supplied by the programmer and may be any 
hexadecimal integer less than 216 — 1. 


Note: In this case, the entry name must be defined by an EXTRN instruc- 
tion to obtain proper linkage. 

This statement defines the save area needed by the subprogram. The xx 
represents the minimum size of the save area required; however, the pro- 
grammer is advised to use a save area of 18 full-words for all subprograms. 
(The minimum save area requirements are given in Tables 2 through 6 
for the mathematical subprograms and in Table 16 for the service sub- 
programs. ) 

Aligns the argument list at a full-word boundary. 

Indicates the first byte of the only argument. 

Contains the address of the argument. 


Contains the address of the first argument. 
Contains the address of the second argument, 


Indicates the first byte of the last argument. 
Contains the address of the last argument. 


normalized floating-point real or complex value. An 
integer answer is stored in general register 0, a real 
answer is stored in floating-point register 0, and a 
complex answer is stored in floating-point registers 
0 and 2. 

Tables 2 through 6 contain additional information for 
using the mathematical subprograms in an assembler 
language program. These tables give the floating-point 
registers that are used by the subprogram and the 
save area required by the subprogram. 


Service Subprograms 

The service subprograms do not use the floating-point 
registers during execution; however, each service sub- 
program requires a save area. The minimum size of 
the save area depends upon the subprogram to be 
used and is given in Table 16. 


Table 16. Assembler Information for the Service Subprograms 


Subprogram Entry 
Name Name(s) 






Save Area 
(Full Words ) 




























IHCFDUMP DUMP 
PDUMP 
THCFDVCH DVCHK 
IHCFEXIT EXIT 
IHCFOVER OVERFL 
IHCFSLIT SLITE 





SLITET 
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Appendix F. 


Sample Storage Printouts 


Appendix F contains a sample printout for each dump 
format that can be specified for the ncFDuMP subpro- 
gram. The printouts are given in the following order: 
hexadecimal, LOGICAL *1, LOGICAL *4, INTEGER *2, 
INTEGER *4, REAL *4, REAL *8, COMPLEX *8, COMPLEX *16, 
and literal (see Figure 3). 


- CALL PDUMP 


OOA3E0 


006CC8 
Q060Fe 


CALL PDUMP 
OO6E1E 
CALL PDUMP 
006E10 
CALL PDUMP 
006E18 
OO6ELA 
006E1C 
CALL PDUMP 


006E20 
006E48 


CALL PDUMP 


QO6eECO 


CALL PDUMP 
006DC8 


CALL PDUMP 


006CD0 
CALL PUUMP 


OO6CEO 
CALL PDUMP 


WITH HEXADEC IMAL FORMAT SPECIFIED 
485F5E10 00000000 485F5E10 10000000 


42B00000 00000000 00000000 00000000 
CQ000000 10000000 41200000 41566666 


WITH LOGICAL*1 FORMAT SPECIFIED 
T OF 
WITH LOGICAL*4 FORMAT SPECIFIED 
FOOT 
wITH INTEGER*2 FORMAT SPECIFIED 
10 
-100 
10 
wITH INTEGER¥4 FORMAT SPECIFIED 


1 2 3 
11 12 


WITH REAL*4 FORMAT SPECIFIED 
0. 20Q000000E ol 0-53999996E O01 


WITH REAL*8 FORMAT SPECIFIED 
0.1759999559599995D 03 
WITH COMPLEX*8 FORMAT SPECIFIED 


(3.0000000+4.0000000) 
WITH COMPLEX*1L6 FORMAT SPECIFIED 


(0.59999999999999990 10. 9999999999999990 } 


WITH LITERAL FORMAT SPECIFIED 


006E50 THIS ARRAY CONTAINS ALPHAMERIC DATA 


Figure 3. Sample Storage Printouts 
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42100000 


00000000 00000000 00000000 00000000 00000000 
0000000C 41100000 


4 5 6 7 


(420000000, &.0000000) 


00000000 00000000 


(-0.9999999999999990 »-0.9999999999999990 ) 


00000000 


10 
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WISE oars cea Shenoy weak Seas So he eee oe ee cena 9 
IHCSERF subprogram 
AMOFNN s, 4.9.1259-08 Ot. 55s ality the docs Meee 31-32 
effect of an argument error............. 0000.0. ee eaee 32 
WCTIOnMANCe . hae cia ede As cca te dns catia: 42 
SIZE rit Ae ae Geis Bre UN Oe cuore iad fi eee de 48 
MISO3 opt Sue Od. ected oat ne een in tage dias eet Mine a area ante 12 
IHCSEXP subprogram 
al@orithn cs oe hose Ge ete ae ty ee ef Ree ics 32 
effect of an argument error.......... 00.00 c eee eee 33 
@ITON: MESSALO: 5 345 RRs Kes weeks yes dna gelato os 45 
DCTIONMNANCE” cl bl Mog h Put oak in oivaad eed Sako aids 42 
GILES 2 Sieh Sorte Web ed: Ria Hoe Rag Gy RN AO OSI er a ae eS 48 
LSC) in, ay hee ak er ay ee et AG a hte ne eerie iene oe 8 
IHCSGAMA subprogram 
IPOVIC EI 2 See a cia at Belg Sls Gh Cake are a tl Mane Ne gt 33 
effect of an argument error.................0.0 000s 33 
CTLOE “MESSAMES fo wete ans Gb Pee eds Se hae bas AT 
PRTIOLMNANCE: iii eh oe hs GENE ee 39, 42 
AZO ck. 3 Ble GSAS alt ase OPEL BRA Alach tes ehere Senco tah area dg 48 
MISE rhs otlcet he Soa ate at eh ene ate Gaeta Hand 29 Giedt 8 12 
IHCSLOG subprogram 
AISOMEDT: 35058 dS oue naa hehe nat a Palen Ae ea 33-34 
effect of an argument error........... 00.00. c eeu 34 
GXTOL INESSATO 5 5h e's Wa ee ERIS ha Ae RRs 45 
PELLOLMMADCEs co n% 6.0 25s oe Soe ew FES RA Rat See 39 
SIZE 5 ok hea ea Ray See ERE OEE Na ee OES 48 
MISO). <2. 5 ok aesaree Se need SS ages oes OLA 8 
IHCSSCN subprogram 
GONE, ake cs ace soe Rae he Bae hes. ee 34-35 
effect of an argument error.............. 00000 e eee 35 
€TTOL MESSATE So cek oF Re a ae | A5 
PETIOPMANCE. of. :h6156 54 WR ee Rs oa dean aad whe eho 40, 42 
BIDE ks sos dt Araceae Rin id ulastioe alterna ie hor a gel whsagt hea Be nee ee 48 
WSO. aie 4 cso SAE A LEAE Eee SERENE ee eee 10 
IHCSSCNH subprogram 
ALGOMITAINY, 6 hor oye Oe a ee hee hele Rares gta thaee hat a 35 
effect of an argument error............. 000000 e eee 35 
CrrOr MESSAGE! Pea eee ete aba aa SR OF 45 
POLFOPMANCE? {4.6.6 pee WS Sha ee VE is kk eG 40, 42 
SEE 3 canna SA Re BAA a ee eh A a iste a ala. eH A Sued oxe eh ge 48 
MSOY Sh a Wet hee hid eel Diy. baa ln, Bab eset Aaa a Il 
IHCSSQRT subprogram 
BI GOPU I e424, 31ers ble a eee. ool te eens 35-36 


effect of an argument error..................0..0.000. 36 


eEYror MESSAGE .... 6... ee ence eens 45 
PETlOrMAanCe.:. 4.50 Sk Be hing pus we eae See etyne Ba ee 42 
S1ZO ish sie AE ech Sen Sok AEE EE, Lik Wak Serge PR WE aS 48 
USE) 35 5 ca eon ee bd 8G OR 8 te Send GHB ad wi oie 8 

IHCSTANH subprogram 
BIQOVIENIN Sues olkets heeyie 2 oeatce ce he ol re See 36 
effect of an argument error..................0.0.0000, 36 
PETLOVIMANCE? 5 skis het hak nk Bede heen au 43 
SIZE: af suds cosets Cand saa nachedaed an Ag Shee oan a ee Kean eae oe es Oke 48 
USO) iis Side ek aad ae Beso SiR Sandth See Mbesebla aoe CANN 11 

IHCSTNCT subprogram 
SIWORILINN, cere ee acerap lard Md Reon See Os ek aS 36-37 
effect of an argument error...............0 00.0000 eae 37 
EYYOr MESSAGES .. 1. ee eee 45 
performance .........0.. 0.000 e eee cence ees 40, 43 

TSIZO i a! Nak Sia wie, Auremegew tk ean phd cone yh OIE ats in Resin 48 
USE eee kaise Lee yee se ore a cae Slims to ae ccaoeae BNE he Mere ed acd a eb ad ca 10 

INT (see IHCFIFIX) 

Implicity called subprograms.................... 5, 6, 14-15 
List Ob-2 op Seo 0d caw hack She ea Se eee nbn oe eee 14 
TOSUICOE USE 2.44 Co scbans 2 he esto ag anne paipeleee 2A 15 
SIZC oioticies es Both yea Ge ait ete Se ER SW ee We NS Se Read a cede 48 
MISC sors aceasta eed aha etary ent Bee asa ar ae ae 14 

Interruption procedures ............0 00.000 e eee eee 44 

Linkage egtor 4 2.c)-coohas 6.2984 RSA S oe KSSH Le REIS A 5 

Logarithmic subprograms........... 7, 8, 19, 25-26, 33-34, 48 

Log-gamma subprograms ........... ...... 7, 12, 25, 33, 48 

Machine indicator test subprograms............... 16, 49, 51 

Mathematical subprograms ................. 00.000 e eee 5, 6 
WLBOFILNINS *y5.ccaceyh di pt Diced 4a dash ioe Ss ne aed eee eae 18-37 
CERNOR 4 2b he ree re oe heath aes eee Peek es 5 
explicitly called 06. ost0¢20 win ded oda Sa ie edd es ed 6-13 
implicitly Called g:5 6:4 cece Bae weeks ER ee SL ORR 14-15 
MISE Ob oir a ashe at atet nd od, B-Series ea ATA 7,14 
DETIOFMANCE. 6.604 2b ead oo eke, Les 34 EG Aaa SOEs 38-43 
SIZES: al iste Se die Be a rea tet ee ee cans, hc eae nea 48 
uisein FORTRAN 234 ooh swe 4 ee ve te Lake pe tet as 6-15 
use in assembler language.................-0-0005. 50-51 

Maximum value subprograms .................. 7, 12-13, 48 
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MAXO (see IHCFMAXI) 
MAXI (see IHCFMAXR) 
MINO (see IHCFMAXI) 
MINI (see IHCFMAXR) 


Minimum value subprograms................... 7, 12-13, 48 
MOD (see IHCFMODI) 

Modular arithmetic subprograms .................. 7, 13, 48 
Natural logarithm subprograms..... . 7, 8, 19, 25-26, 33-34, 48 


OVERFL (see IHCFOVER) 
PDUMP (see IHCFDUMP) 


Psendo sense lights sss ie it katana beh ee SE Sale cedars 16 
Relative error: << cot uae ns Seta Wea hehe be arches 18, 38 
Sample dump printouts ............0 00.000 c cece ees 52 
Sampling techniques ..............0 0000 cece cence eee 38 
Sense: lights: 4.52.4 aaa Paaa aaa eR ars hee Chae 16-17 
Service subprograms 
machine indicator test ...........0 0.000000 c ccc 16 
BAZ ib tiatns ht de SHS ke gh satin GYoocat. & dnten Op ance tog gare cn a ».. 49 
use in assembler language...................00000- 50-51 
use in FORTRAN .............. 0.0: e eee eee ees 16-17 
RUEUEE Ye ics oc beeen s ok Be ws AR Taf Bs aca pepe land area 16-17 
SIN (see IHCSSCN ) 
Sine subprograms............ 7, 9-10, 20-21, 26-27, 34-35, 48 


SINH (see IHCSSCNH ) 
SLITE (see IHCFSLIT) 
SLITET (see IHCFSLIT) 
SQRT (see IHCSSQRT ) 


Square root subprograms.......... 7, 8, 19-20, 27, 35-36, 48 
Standard deviation ... 0.0.00... 00.0: cece eee nee 38 
Storage estimates ... 0.0... ccc ee eens 48-49 
Storage printouts .........0. 0.0, c cece cece ee eens 52 
TAN (see IHCSTNCT) 

Tangent subprograms................ 7, 10, 28-29, 36-37, 48 
TANH (see IHCSTANH ) 

Timing statistics ...... 0.000. .c ccc eee 38-48 
Trigonometric subprograms . .7, 9-10, 20-23, 26-31, 34-37, 48 
Truncation subprograms ...............0 00000 eae 7, 13, 48 
Utility subprograms .................. 0045. 16-17, 48, 50-51 
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